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An explicitly orbital-dependent correlation energy functional is proposed, which is to be used in 
combination with the orbital-dependent exchange energy functional in energy-band calculations. It 
bears a close resemblance to the second-order direct and exchange perturbation terms calculated 
with Kohn-Sham orbitals and Kohn-Sham energies except that one of the two Coulomb interactions 
entering each term is replaced by an effective interaction which contains information about long-, 
intermediate-, and short-range correlations beyond second-order perturbation theory. Such an ef- 
fective interaction can rigorously be defined for the correlation energy of the uniform electron liquid 
and is evaluated with high accuracy in order to apply to the orbital-dependent correlation energy 
functional. The coupling-constant-averaged spin-parallel and spin-antiparallel pair correlation func- 
tions are also evaluated with high accuracy for the electron liquid. The present orbital-dependent 
correlation energy functional with the effective interaction borrowed from the electron liquid is 
valid for tightly-binding electrons as well as for nearly-free electrons in marked contrast with the 
conventional local density approximation. 
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I. INTRODUCTION 

The band theory of solids has been developed with the 
help of the knowledge of exchange and correlation ac- 
quired from the study of the uniform electron liquid. The 
exchange difficulty arising from Hartree-Fock equations 
has been overcome by two different methods. One is the 
standard many-body theory which has established the 
quasi-particle picture through Dyson equations involv- 
ing the self-energy operator responsible for all exchange 
and correlation effects on quasi-particles and clarified the 
importance of dynamically screening the exchange po- 
tential in the evaluation of the quasi-particle energy dis- 
persion of the electron liquid Q, The other is the 
Slater's intuition 0, Q that the pathological behavior 
of the wavenumber-dependence of the Hartree-Fock ex- 
change potential for the electron liquid around the Fermi 
wavenumber should be smoothed out by averaging the 
potential over all occupied states and that real many- 
electron systems may be regarded as a locally uniform 
electron liquid with the same local electron density to 
employ this averaged exchange potential. The Slater's 
intuition has developed into a rigorously founded many- 
body theory valid for the ground state, namely, density- 
functional theory (DFT) HQ. 

In this paper we propose a new method for further de- 
veloping the band theory along DFT. For this purpose 
we borrow detailed knowledge of long-, intermediate- 
, and short-range correlations from the electron liquid, 
not such averaged knowledge as the magnitude of its 
exchange-correlation energy on which the conventional 
local density approximation (LDA) is based. First we 
give a new expression for the correlation energy of the 
electron liquid since it is needed for a proposal of an 



orbital-dependent correlation energy functional for use 
in energy-band calculations. It resembles second-order 
direct and exchange perturbation terms, but one of the 
two Coulomb interactions entering each term is replaced 
by an effective interaction which contains all effects from 
the higher-order perturbation terms beyond second-order 
perturbation theory. The effective interaction thus de- 
fined is reduced in magnitude from the bare Coulomb 
interaction for all distances since it takes full account 
of long-, intermediate-, and short-range correlations in 
agreement with the Pauli principle. We evaluate the ef- 
fective interaction with high accuracy over the entire re- 
gion of metallic densities. 

According to Moruzzi, Janak, and Williams Q, the 
lowest possible uniform density realized in the interstitial 
region outside muffin-tin spheres of metals corresponds 
approximately to the critical density (r s — 5.25) where 
the compressibility of the electron liquid becomes diver- 
gent, indicative of an instability of this model. In the 
construction of an orbital-dependent correlation energy 
functional for energy-band calculations, we think it is 
physically sound to borrow the knowledge of the effec- 
tive interaction from the electron liquid above the lowest 
critical density where the system remains thermodynami- 
cally stable. Note that the LDA regards real systems as a 
locally uniform electron liquid and borrows its exchange- 
correlation energy even below the critical density. 

The orbital-dependent correlation energy functional we 
propose here consists of a direct and exchange pair of 
second-order perturbation like terms constructed with 
Kohn-Sham orbitals and Kohn-Sham energies, in which 
one of the two Coulomb interactions in each term is re- 
placed by the effective interaction borrowed from the 
electron liquid. A combination of this orbital-dependent 
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correlation energy functional and the orbital-dependent 
exchange energy functional is of course reduced to the 
LDA in the limit of uniform density. The use of the 
orbital-dependent exchange energy functional in DFT 
makes it possible to cancel exactly the spurious self- 
interaction terms involved in the classical Hartree energy 
functional. Furthermore, the present orbital-dependent 
correlation energy functional cancels exactly higher-order 
self-interaction terms since it consists of a direct and ex- 
change pair of terms. We would like to stress that the 
present orbital-dependent correlation energy functional 
reflects the anti-symmetric property of the many-electron 
wave function in a similar way as the Hartree and the 
exchange energy functionals in the Hartree-Fock approx- 
imation. 

In order to obtain the exchange and correlation po- 
tential v xc (r) in DFT, it is necessary to evaluate the 
functional derivative of the exchange and correlation en- 
ergy functional E xc with respect to the electron density 
n(r). However, the present correlation energy functional 
as well as the exact exchange energy functional is an ex- 
plicit functional of Kohn-Sham orbitals, and hence is an 
implicit functional of n(r). The method of evaluating 
the exchange potential through the functional derivative 
of the orbital-dependent exchange energy functional with 
respect to Kohn-Sham orbitals has been established and 
termed the optimized effective potential method (OEP) 

The present correlation energy functional is given as 
a functional of both Kohn-Sham orbitals and the effec- 
tive interaction borrowed from the electron liquid. Then, 
the OEP method can be applied to the evaluation of the 
correlation potential but the functional derivative of the 
effective interaction with respect to the electron density 
n(r) has to be treated separately by resorting to an ap- 
propriate approximation. 

Another problem is how to choose the density parame- 
ter for the effective potential we borrow from the electron 
liquid. In this respect we may utilize the fact that the 
compressibility of the electron liquid calculated for the 
uniform density in the interstitial region outside muffin- 
tin spheres is generally a good approximation to the com- 
pressibility of real metals. Hence it is a reasonable ap- 
proximation to choose as the density parameter the value 
appropriate for the uniform electron density in the inter- 
stitial region of metals calculated from the LDA. The 
best choice of the density parameter is the optimization 
that minimizes the calculated ground state energy in the 
present theory as a function of the density parameter. 

One of the remarkable merits of DFT is that the 
correlation-induced reconstruction of Kohn-Sham or- 
bitals accompanied by a change in the electron density 
n(r) of the system can in principle be described with the 
local correlation potential v c (r). This can be realized 
only when accurate information about orbital-dependent 
correlation energy functional £ c [{^?j}] and the poten- 
tial v c (r) defined by 5E c [{<pi}]/5n(r) is avaiable in ad- 
dition to the orbital-dependent exchange energy func- 



tional E x [{ifi}] and the resulting exact exchange poten- 
tial v x (r). 

This reconstruction of Kohn-Sham orbitals is beyond 
the scope of the LDA since the LDA regards real sys- 
tems as a locally uniform electron liquid and borrows 
its density dependence of the exchange-correlation en- 
ergy. The presence of spurious self-interaction terms and 
the absence of the correlation-induced reconstruction of 
Kohn-Sham orbitals are the main two reasons why the 
LDA fails to give the correct evaluation of the band gap 
of semiconductors and the energy-band structure of the 
so-called strongly correlated electron systems. 

In the present theory, the exchange and correlation 
energy functionals are both orbital-dependent and at the 
same time they are free from self-interaction errors. The 
present correlation potential v c (r) depends implicitly on 
the electron density n(r') at the neighboring position r'of 
r through the orbital-dependence of £^[{1^}]. 

An important prediction the present theory makes is 
that correlation may reduce the screening of the nuclei 
at short distances to produce a shrinkage of the elec- 
tron distribution around individual nuclei of the system 
and in accordance with this redistribution of the electron 
density Kohn-Sham orbitals may be reconstructed par- 
ticularly in the neighborhood of the Fermi level or the 
energy gap. 

An accurate description of exchange-correlation ef- 
fects on the energy-band structure is expected from the 
present theory. Its orbital-dependent correlation energy 
functional, if applied to the so-called strongly correlated 
electron systems, will give rise to significant differences 
in the energy-band structure compared with the case of 
the orbital-dependent exchange energy functional alone, 
particularly in the neighborhood of the Fermi level or the 
energy gap. 

Since the early stage of the development in DFT a 
number of authors have been interested in the rela- 
tionship between Kohn-Sham equations and Dyson equa- 
tions. DFT is an exact many-body theory for the ground 
state in the framework of the self-consistent one-electron 
theory. The Hohenberg-Kohn variational principle leads 
to a set of Kohn-Sham equations describing the reference 
non-interacting system with the same electron density 
as the real many-electron system. In this reference sys- 
tem electrons interact only through the local exchange- 
correlation potential v xc (r). On the other hand, the stan- 
dard many-body theory is founded on the variational 
principle [l(J, uM t na t the ground state energy defined 
as a functional of the many-body one-electron Green's 
function G(r, r', oj) is stationary with respect to any vari- 
ation in G(r, r', lj). This stationary property leads to 
Dyson equations for the one-electron Green's function. 
The one-electron Green's function has the well-founded 
physical meaning to give the concept of quasi-particles 
as low-lying elementary excitations of extended many- 
electron systems. 

To make a comparison between Dyson equations and 
Kohn-Sham equations we have evaluated with high accu- 



racy the quasi-particle energy dispersion of the electron 
liquid over the entire region of metallic densities. It is 
found that in the limit of uniform density energy spectra 
of Kohn-Sham equations are a rather good approxima- 
tion to those of Dyson equations unless they exceed a 
critical energy beyond which screening no longer works 
around the Fermi energy plus the plasmon excitation en- 
ergy. 

In order to justify the application of Kohn-Sham equa- 
tions to the band theory it is pointed out that the implicit 
dependence of the self-energy £(r, r', oj\ [G]) in Dyson 
equations on the self-consistent Green's function G is in 
one-to-one correspondence to the implicit dependence of 
exchange-correlation potential v xc (r; [n]) in Kohn-Sham 
equations on the self-consistent electron density n(r). 

The present paper is organized as follows. In Sec. II 
we give a detailed analysis of the dielectric formulation 
for the study of the electron liquid. This is a section 
requisite for the introduction of a new expression for the 
correlation energy. In Sec. Ill we calculate the effec- 
tive interaction by interpolating between long-range cor- 
relation in the random phase approximation (RPA) and 
short-range correlation in the particle-particle ladder ap- 
proximation such that the exchange counterparts are self- 
consistently included. We also calculate the coupling- 
constant-averaged spin-parallel and spin-antiparallel pair 
correlation functions as well as the correlation energy. In 
Sec. IV we propose an orbital-dependent correlation en- 
ergy functional for use in energy-band calculations. Some 
characteristic properties of this correlation energy func- 
tional and how to use it in practical calculations are dis- 
cussed. In Sec. V we show an accurate calculation of the 
quasi-particle energy dispersion of the electron liquid. In 
the last section we give some concluding remarks. 



II. THEORY OF ELECTRON LIQUID 

A. Exchange correction to the RPA 

The random phase approximation (RPA) 0, 0, IbH 
Il5j is a good starting point to describe long-range corre- 
lation characteristic of the electron liquid, but it violates 
the Pauli principle since the second- and higher-order 
exchange counterparts of the RPA series of perturba- 
tion terms are missing. Hubbard [l4l fl5| was the first 
who made an attempt to allow for these exchange cor- 
rections. His approximation is an interpolation between 
long-range correlation in the RPA and short-range cor- 
relation in second-order perturbation theory, but all the 
coefficients of third- and higher-order exchange terms in 
his approximation are miscounted. The reason for this 
miscounting is that he failed to enumerate all the pos- 
sible proper polarization diagrams leading to the same 
third- and higher-order exchange terms. This is not seri- 
ous for his evaluation of the ground state energy E g , but 
is very serious for the purpose of evaluating the quasipar- 
ticle energy dispersion E(p) and the Landau interaction 




FIG. 1: Second- and third-order Goldstone energy diagrams; 
upward and downward solid lines denote particles and holes, 
respectively, and dotted lines denote Coulomb interactions. 
(2a) and (2b): a direct and exchange pair of second-order di- 
agrams; (3a): third-order PRA diagram; (3b), (3c) and (3d): 
singly exchanged diagrams; (3e),(3f) and (3g): doubly ex- 
changed diagrams; (3h): triply exchanged diagram; (3i) and 
(3j): third-order diagrams with a Coulomb interaction line 
inserted between two particle or two hole lines in the second- 
order direct diagram; (3k) and (31) : exchange counterparts of 
(3i) and (3j). 



function f aa (p,p') which are closely related to the first 
and the second functional derivations of the ground state 
energy E g [G] with respect to the quasi-particle distribu- 
tion function n(p), respectively fig . 

By way of illustration we consider the exchange coun- 
terparts of the third-order PRA diagram (see Fig. 1). 
Three different exchange counterparts can be obtained by 
exchanging two incoming (or outgoing) particle or hole 
lines attached to one of the three Coulomb interaction 
lines in the RPA diagram. Additional three different 
exchange counterparts can also be obtained by doubly 
exchanging the PRA diagram. The last exchange coun- 
terpart can be obtained by exchanging all of the three 
interaction lines; this is a diagram with one interaction 
line inserted between a pair of particle and hole belonging 
to different bubbles in the second-order direct diagram. 
The Hubbard approximation miscounts the coefficients 
of singly and doubly exchanged perturbation terms and 
furthermore neglects the triply exchanged perturbation 
term. Besides the eight third-order diagrams we have just 
enumerated, there remain more vertex-type diagrams of 
the same order which are not derivable from the same 
order RPA diagram by exchanging, i.e., two diagrams 
with one Coulomb interaction line inserted between two 
particles or two holes in the second-order direct diagram 
and their exchange counterparts. Finally, there are third- 
order diagrams with a first-order self-energy correction 
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inserted into the second-order direct and exchange dia- 
grams. 

Even if all the possible exchange counterparts of the 
RPA series are correctly calculated, another problem 
arises instead. The resulting approximation is unable 
to give the correct description of long-range correlation 
characteristic of the electron liquid. 

The conserving approximation method by Baym and 
Kadanoff 0,0] can raise levels of approximations main- 
taining the correct description of long-range correlation 
by improving the generating functional $ progressively. 
This method guarantees conservation of the local elec- 
tron number, the local electron momentum and the local 
electron energy as well as the f-sum rule. However, it 
does not fulfill the Pauli principle unless the theory is 
exact. As a consequence the requirement of the Pauli 
principle in the framework of the conserving approxima- 
tion method can be a guiding principle leading to the 
exact many-electron theory. 



B. Analysis of the dielectric formulation 



The dielectric formulation is originally suitable for the 
description of long-range correlation of the electron liq- 
uid. It discriminates between proper and improper polar- 
ization functions. Improper polarization functions are re- 
sponsible for long-range correlation and proper polariza- 
tion functions for short-range correlation and exchange. 
By a diagrammatic analysis we shall investigate how the 
requirement of the Pauli principle should be fulfilled in 
the dielectric formulation and thereby give the integral 
equation for the particle-hole irreducible interaction in a 
diagrammatic representation (see Fig. 2). 

We shall start with the following equations 0, 120I l2ll 
I22I |23| for the spin-parallel and spin-antiparallel parts of 
the static structure factor 5(q). 



= ^E(^=i^ <^ " Im {^)( l -^)) + '• (<l '" , 

irnj { 2 2 e(q,w) J 

h r° fl l ^(q,^(q)7r(q,^) 

= — / dulml -7r 6 (q,w) - 

nn J [2 2 e(q, uj) 



(1) 



(2) 



where (■ • ■ ) denotes the expectation value with respect 
to the exact ground state and is the density fluctu- 
ation operator; the dielectric function e(q, uj) is defined 
as e(q, uj) = 1 + u(q)7r(q, uj) and w(q) = 47re 2 /q 2 . The 
proper polarization function 7r(q, uj) is divided into ex- 
change and direct parts as 7r(q, uj) = 7r a (q, uj) + 7T(,(q,a;) 
(see Fig. 2 (a), (b) and (c)). The static structure factor 
5(q) that gives the description of charge-fluctuations is 
written as 



5(q) - S TT (q)+S U (q) 

= ^<(pi+^)(P T - q +^- q )> 

= AfW 1 U 1 )). (3 ) 



On the other hand, the exchange part of S(q), S(q) that 



gives the description of spin-fluctuations is written as 
5(q) - STT (q) _ 5 U( q) 

= ^((Pq-P q )(P- q -^q)) 

= — / dui Im{w a (q, uj)}. (4) 

The exchange function 7r a (q, uj) comes from those di- 
agrams in which two external points are connected 
through particle and hole lines; the function 7r a (q, uj) 
is spin-dependent and makes a direct contribution to 
S'TT(q) alone; it is reduced to the Lindhard function 
7ro(q, uj) in the non-interacting limit. The direct func- 
tion 7Tb(q, uj), on the other hand, comes from those di- 
agrams in which two external points are not connected 
through particle and hole lines; the function 7T£,(q, to) is 
spin-independent and makes a contribution to S'' (q) and 
S'^(q) equivalently. The improper polarization function, 
7r(q, w)w(q)7r(q, ui)/e(q, uj), whose diagrams can be sepa- 
rated into two pieces by cutting a single Coulomb inter- 
action line in them, also makes a contribution to S'^(q) 
and 5^1 (q) equivalently. 
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FIG. 2: A diagrammatic representation of the integral equa- 
tion for the particle- hole irreducible interaction, (a): the in- 
tegral equation for the screened interaction; wavy and dot- 
ted lines denote the screened and the bare Coulomb inter- 
actions, respectively, and a bubble with a cross-hatched tri- 
angle denotes the proper polarization function, (b): the 
integral equation for the proper vertex part; an ellipse de- 
notes the particle-hole irreducible interaction, (c): decom- 
position of the particle-hole irreducible interaction into two 
parts, direct and exchange, (d): those polarization diagrams 
which integrate to zero due to the Pauli principle; the first 
infinite series of diagrams denotes the exchange part of the 
proper polarization function excluding the Lindhard function, 
7T a (<7, to) — 7ro(<7, to); the second infinite series of diagrams with 
a factor of 1 /2 denotes the direct part of the proper polariza- 
tion function, Hb(q,u>) and the third part, two bubbles con- 
nected with a wavy line denotes the improper polarization 
function, (e) and (f): the integral equations for the exchange 
and direct parts of the particle-hole irreducible interaction, 
respectively, (g): the integral equation for the total particle- 
hole irreducible interaction; the right-hand side consists of 
three parts, improper, doubly irreducible, and proper and re- 
ducible, (h): decomposition of the doubly irreducible inter- 
action into two parts, exchange and direct; a square denotes 
the doubly irreducible interaction. 



It is important to note that the exchange function 
7r a (q, uj) makes a contribution to S'T'r(q) and S'^(q) indi- 
rectly through the improper polarization function com- 
mon to the two. In other words, the two functions 5^(q) 
and .S^-l^q) are coupled through 7r Q (q, uj) in the improper 
polarization function. This means that there is correla- 
tion between charge-fluctuations and spin-fluctuations. 

The proper polarization function 7r(q, to) that is the 
key quantity in the dielectric formulation is given as 

7r(q,o;) = -±£ / ^(p, e)G a (p + q, e + u) 
x A(p,e;q,w), 5^+0, (5) 



where G cr (p, e) is the one-electron Green's function and 
the proper vertex part A(p, e; q, uj) satisfies the following 
integral equation. 



A(p,e;q,w) 

n 



de' 



x G a , (p', e')G a < (p' + q, e' + u>)A(p', e'; q, «). (6) 

The particle-hole irreducible interaction / <T<T (p, e;p',e'; 
q, uj) can be divided into two parts, direct and ex- 
change (see Fig. 2 (c)). The function 7r a (q, uj), apart 
from 7To(q, to), is constructed only from exchange parts, 
whereas the function 7Tf, (q, uj) is constructed from both 
direct and exchange parts but must include at least one 
direct part. 

The spin-parallel and spin-antiparallel pair correlation 
functions are given respectively as 
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n 



to = 1 + ^E^ TT (q)-i) 



to = i + ^E^ u (qK q - r - 



(7) 

(8) 



The Pauli principle means that g'"(0) = 0. This is sat- 
isfied only by 7To(q, u>) or its modification |l9j| in which 
two non-interacting one-electron Green's functions are re- 
placed by two exact one-electron Green's functions. The 
Pauli principle then means that 



f°° ( 1 

V / dujlml 7r Q (q,w) - 7r (q,w) + -7r b (q, 

q J( > 1 1 

1 7r(q,w)w(q)7r(q,o;) 



e(q,w) 



= 0. 



(9) 



We shall see how a direct term and its exchange coun- 
terpart which exactly cancel each other by integrat- 
ing over q and uj appear in pairs in the equation 
above. First of all, the quantity n a (q,uj) — 7r (q, uj) is 
nothing but the exchange counterpart of ^{^(q, to) — 
7r(q,w)'u(q)7r(q,w)/e(q,a;)} (see Fig. 2 (d)). Con- 
sider in detail the interrelation between direct and ex- 
change terms. For this purpose we shall define the 
order of 7r a (q, to) and 7r&(q, to) according to the num- 
ber of irreducible particle-hole interactions included in 
them (see Fig. 2 (a), (b) and (c)). It can easily be 
seen that the exchange counterpart of the improper term 
— ^7r(q, w)u(q)7r(q, w)/e(q, oj) is part of the first-order 
contribution of 7r a (q, to) that consists of a single exchange 
part of I C " T . The exchange counterpart of ^^(q, ui) is 
either the first- or higher-order contributions of 7r a (q, uj). 
Consider the first-order contribution of ^7rfc(q, uj) that 

has a single direct part of l aa . The interaction in it 
can be separated into two parts: One is the doubly irre- 
ducible interaction and the other is the direct mate of all 
the second- and higher-order contributions of 7r a (q, uj) 
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(see Fig. 2 (f)). By doubly irreducible we mean that 
the diagrams have no particle-hole pair in the interme- 
diate states when they are seen not only from bottom 
to top, but also from left to right. All the second- and 
higher-order contributions of s7Tf,(q, uS) include at least 
one direct part of I a(T . The exchange counterpart of all 
these contributions is part of the first-order contribution 
of vr a (q, uj). 

We have seen how various parts of polarization func- 
tions, when integrated over uj and q, should be cancelled 
to fulfill the Pauli principle. From the analysis above we 
may conclude that the exchange part of I aa consists of 
(i) the exchange counterpart of the interaction included 
in the improper polarization function, (ii) the exchange 
part of the doubly irreducible interaction, and (iii) the 
exchange counterpart of the interaction included in all 
the second- and higher-order contributions of ^7T(,(q, uj) 
(see Fig. 2 (e)). We may also conclude that the di- 
rect part of I acr consists of (i) the direct part of the 
doubly irreducible interaction and (ii) the direct mate of 
the interaction included in all the second- and higher- 
order contributions of 7r a (q, uj) (see Fig. 2 (f)). We have 
thus arrived at the same conclusion concerning the func- 
tional structure of the particle-hole irreducible interac- 
tion that Yasuhara and Takada |Q| obtained from the 
exact analysis of the self-energy £o-(p, e) as a functional 
of G CT ,(p',e'). 

Screening, or long-range correlation in the electron liq- 
uid is ascribed to the long-range nature of the Coulomb 
interaction between electrons. An infinite summation 
of many-body perturbation terms is needed for the de- 
scription of screening. As a result, the requirement of 
the Pauli principle in the dielectric formulation is ful- 
filled by the very complicated form of integral equation 
for the particle-hole irreducible interaction diagrammati- 
cally shown in Fig. 2 (g). Thus the unified description of 
screening and exchange is very much difficult to achieve 
in the electron liquid. 

The most important part of the doubly irreducible in- 
teraction shown in Fig. 2 (h) is considered to be an 
infinite series of particle-particle and hole-hole ladder 
interactions (direct and exchange) with respect to the 
screened Coulomb interaction. The inclusion of an infi- 
nite series of particle-particle ladder interactions via the 
bare Coulomb interaction in the many-body perturba- 
tion expansion is indispensable to the proper description 
of short-range correlation as well as to the fulfillment of 
the cusp condition. 

C. Particle- hole irreducible interaction 

The integral equation for l aa (p, e; p', e'; q, uj) we have 
just obtained from a diagrammatic analysis of polariza- 
tion functions (see Fig. 2 (g)) is quite the same as the in- 
tegral equation that Yasuhara and Takada [24| obtained 
from a diagrammatic analysis of the self-energy S CT (p, e) 
as a functional of G a >(p',e') in their derivation of the 



symmetric expression for £ CT (p, e) in which each con- 
stituent G a > of £ CT (p, e) enters in an equivalent way (see 
Ea. (|56|) ). According to many-body theory, this coinci- 
dence is natural since the particle-hole irreducible inter- 
action l aa (p, e; p', e'; q, uj) in the limit q = uj = is iden- 
tical to the functional derivative <5£ CT (p, e)/#G CT '(p', e'). 
Hence the functional structure of I c " 7 (p, e; p', e'; q, uj) is 
identical to that of 5Y lcr (p,e)/8G r ?>(p' ,e'). 

The conserving approximation method by Baym and 
Kadanoff 0,0] requires that an expression for £er(p, e) 
and the one for I aa (p, e; p', e'; q, uj) to be used in the con- 
struction of an approximate proper polarization function 
7r(q, uj) be both ^-derivable. / <T<T (p, e; p', e'; q, uj) has the 
same functional structure as S CT (p, e)/G (J ' (p', e'). The 
conserving approximation method then fulfills one of the 
rigorous self-consistent requirements imposed on the ex- 
act many-body theory. 

With the exact generating functional <£> one can ob- 
tain the exact self-energy and the exact particle-hole ir- 
reducible interaction from its first and second functional 
derivatives, respectively. The standard expression for the 
self-energy is given as 

Mp,e) = -^E / |^G CT (p + q,e + a,) 

x -p^) A( P> e;q,w), 6^+0. (10) 
e(q,w) 

and the one-electron Green's function G CT (p, e) is related 
to the self-energy as 

[G a (p, e)]- 1 = [G° (p, e)]- 1 - E ff (p, e), (11) 

where GjJ.(p, e) is the non-interacting one-electron 
Green's function. Note that ?r(q, ui) and A(p, e;q, us) 
in Eq. (|10|) satisfy Eqs. JSJ and ©, respectively. We 
would like to stress here the importance of satisfying the 
integral equation for the particle-hole irreducible inter- 
action I' 7 ' 7 (p, e; p', e'; q, u>) in the self-consistent evalu- 
ation of various quantities such as G CT (p, e), S a (p,e), 
A(p, e;q, us), 7r(q, uj) and e(q, uj). The accurate knowl- 
edge of I aa (p, e; p', e'; q, uj) makes it possible to evaluate 
all these quantities self-consistently in agreement with 
the Pauli principle. In this sense the integral equation 
for the particle-hole irreducible interaction is the crux of 
many-body theory. 

We have thus clarified that the requirement of the Pauli 
principle in the framework of the dielectric formulation 
can lead to the exact integral equation for the particle- 
hole irreducible interaction I aa (p, e; p', e'; q, uj) and that 
the same requirement imposed on the conserving approx- 
imation method can also lead to the exact many-electron 
theory. 

As can be seen in standard text books (2^, Hf| |23, H3| 

on many-body theory, there are two different expressions 
for the ground state energy: One is written in terms of 
the one-particle Green's function G a (p,e) and the self- 
energy E CT (p, e) and the other in terms of the dynami- 
cally screened Coulomb interaction w(q)/e(q, uj) and the 
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proper polarization function 7r(q, uj); the former can eas- 
ily be transformed into the latter by definition of these 
quantities. The two expressions, if inspected into the 
innermost structure, can be represented with the same 
integral equation for the particle-hole irreducible inter- 
action I aa (p, e; p', e'; q, uj). According to the symmet- 
ric expression for the self-energy E CT (p, e) by Yasuhara 
and Takada [24|, the knowledge of l aa ' (p, e; p', e'; q, uj) 
at q = uj = is sufficient for the evaluation of the ground 
state energy as well as the self-energy Scr(p, e). 

It is important to recognize that the expansion of the 
self-energy S .(p, e) or the proper polarization function 
7r(q, uj) with respect to the screened Coulomb interac- 
tion converges very slowly since correlation is much com- 
plicated in the region of realistic metallic densities. It 
may then be naturally expected that the systematic in- 
clusion of higher-order perturbation terms with respect 
to the screened Coulomb interaction is indispensable for 
the description of strongly correlated features of the elec- 
tron liquid. In fact, there is a great deal of cancella- 
tion due to the Pauli principle occurring among various 
higher-order perturbation terms. For the evaluation of 
the quasi-particle energy dispersion E(p) as well as the 
effective mass m*, it is necessary to consider the inte- 
gral equation for the particle-hole irreducible interaction 
which makes it possible to take systematically into ac- 
count the cancellation between higher-order terms due 
to the Pauli principle. 

Finally we give a brief survey of the present situation 
in the study of the electron liquid. (1) The numerical 
values of the correlation energy calculated from different 
theoretical methods [llllllaillSlllllllllllllllIll 
as well as from the Green's function Monte Carlo method 
HH H0| agree within an accuracy of 0.5 mRy. per elec- 
tron throughout the whole region of metallic densities. 
The compressibility and the spin-susceptibility are also 
evaluated with high accuracy. The RPA overestimates 
the correlation energy in magnitude by about 47% at a 
typical metallic density (r s ~ 4.0) since it lacks higher- 
order particle-particle ladder interactions responsible for 
short-range correlation as well as exchange corrections. 
(2) Accurate evaluations of the spin-parallel and the 
spin-antiparallel pair correlation functions are performed 
[37L l38l l39| . (3) Physical quantities of quasi-particles 
such as the effective mass m* 0, 0, E3- the renor- 
malization factor z |42j , the Landau interaction function 
/ <T<T (p, p' ) [lp a nd the quasiparticle energy dispersion 
E(p) EE 144 l45l (see Sec. V) are all evaluated with 
sufficient accuracy. (4) A self-consistent theory in which 
the integral equation for the particle-hole irreducible in- 
teraction can be fulfilled by iteration in the framework of 
the conserving approximation scheme has been presented 
by Takada [42| and successfully applied to the quanti- 
tative evaluation of the one-electron spectral function 
-Mp, e) 0| and the dynamical structure factor 5(q, uj) 
[49j of the electron liquid at metallic densities. 



III. NEW EXPRESSION FOR CORRELATION 
ENERGY 

A. Definition 

In this section we define a new expression for the cor- 
relation energy of the electron liquid as follows: 

v / q p,(T,p',(T' 



x v(q) 



{l-/(p + q)}{l-/(p'-q)} 



€ p e p+q 



c p -q 



x K//(q)-<Wv e //(-p + p'-q)}, (12) 



where /(p) is the Fermi distrubution function at OK 

h 2 2 

and e p — ~2^-- The expression bears second-order di- 
rect and exchange perturbation terms, but one of the 
two Coulomb interactions in each term is replaced by 
an effective interaction v e ff(q) which contains informa- 
tion about third- and higher-order contributions in the 
perturbation expansion. In Sec. Ill B, C and D we eval- 
uate the effective interaction v e f / (q) by interpolating be- 
tween long-range correlation in the RPA and short-range 
correlation in the particle-particle ladder approximation 
such that the exchange counterpart u e //(— p + p' — q) 
and its feedback effect on the direct interaction v e ff(q) 
are evaluated in a self-consistent way. This expression 
for the correlation energy may be understood from intu- 
ition or from a diagrammatic consideration. Strictly, the 
effective interaction u e //(q) above should be defined as 
v e ff(p, p'; q). Here we employ its averaged value over p 
and p' within the Fermi sphere. 

First we shall see how the above expression for E c 
is related to the dielectric formulation we have de- 
scribed in Sec. II B. The direct term of E c comes 
from the direct part of the total polarization function 
7Th(q, uj) — 7r(q, cj)w(q)7r(q, cj)/e(q, uj) and the exchange 
term of E c from the exchange part of the total polariza- 
tion function excluding the zeroth order Lindhard func- 
tion 7T a (q, uj) - 7T (q, uj). 

We start with the standard expression for the 
exchange-correlation energy, i.e., the integral of the po- 
tential energy over the coupling contant. 



^ = ?£^{S(q)-l}, 



(13) 



where n is the electron density and 5(q) denotes the 
coupling-constant-averaged static structure factor de- 
fined as 



5(q) 



1 



1 



d(e 2 )5(q), S(q) = -< Pq p_ q >. (14) 



The integral over the coupling constant necessarily re- 
duces the magnitude of the potential energy since it im- 
plies the addition of an increase in the kinetic energy 
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caused by correlation to the potential energy. This corre- 
lational increase in the kinetic energy cannot be neglected 
since it cancels about one third of the correlational low- 
ering in the potential energy at a typical metallic density, 
r s = 4.0; in the limit r s — > the correlational lowering 
in the potential energy is by half cancelled by the corre- 
lational increase in the kinetic energy. 

The function S(q) can be decomposed into two parts as 
5(q) = S HF (q) + S c (q). The Hartree-Fock contribution 
S HF (q) is independent of the coupling constant. Then 
the correlation energy is given as 



(15) 



correlation function is defined as 



gH (r) = 1 + ^ { 5TT (q) _ l}( 



(20) 



Since S HF (q) satisfies that g''(0) — 0, the relation 
that X^q^c^q) — has to be satisfied. This identity 
can be easily checked by a simple transformation of the 
wavenumber variables in Eq. H19(l . The Pauli principle 
is automatically fulfilled in the present theory. Instead, 
we have to deal with much difficulty in the self-consistent 
determination of the effective interaction v e f / (q) , as will 
be seen in Sec. Ill B, C and D. 



where S c (q) can be decomposed into direct and exchange 
parts, 5(q) = 5 c d (q) + 5r(q). 

Let us compare Eq. I|15|l with the new expression for 
E c we have defined by Eq. (|12() . Then the functions 
5^ (q) and S^ c (q) may be identified as 



3fo>=nUj EE /(P)/(P') 

P,p' <7,<T' 

„ {l-/(p + q)}{l-/(p'-q)} 



£p ^p+q > ^ 

2 



«e//(q), (16) 



p'-q 



^) = 4uJ EE /(p)/(p') 

p,p' a,cr ; 

„ {l-/(p + q)}{l-/(p'-q)} 



p p+q 



c p -q 



x <W«e//(-P + p' q) 



(17) 



The spin-parallel and spin-antiparallel components of 
5 c (q) are defined as 5^(q)/2 and ^(q)/2 + S° x (q), re- 
spectively, and given as 



£ c u (q) 



i (l 



££/(p)/(p') 



p.p' cr,cr / 

v {l-/(p + q)}{l-/(p'-q)} 

Cp ^p+q "r" ^p' — ^p'— q 

«r«.) = ^(i) 2 EE«p)/(p') 

p.p' <r,a' 

„ {l-/(p + q)}{l-/(p'-q)} 



«e//(q), (18) 



c p+q 



e p' e p'-q 

£ /(p)/(p') 



«e//(q) 



1(1 

P,(7,p' ,£X' 

v {l-/(p + q)}{l-/(p'-q)} 

X 1 

x <W« e //(-p + p' - q)- 



(19) 



Note that this expression for S 1 ^ (q) fulfills the Pauli prin- 
ciple. The coupling-constant-averaged spin-parallel pair 



B. New expression for the RPA correlation energy 

In this subsection we calculate the effective interaction 
v e ff(q) appropriate for the RPA. The correlation contri- 
bution of the static structure factor in the RPA is given 

as 



RPA 



(q) = -S"*(q) 



duAxa- 



7ro(q,w) 



1 + w(q)7r (q,w) 



(21) 



where 7To(q, w) is the Lindhard function. The coupling 
constant integral of S pp (q) can be easily performed as 



HF i 



A f°° fa 1 tan -i ^(q)Im[7r (q,q;)] 
™./o w (q) 1 + w(q)Re[7r (q,cj)]' 



(22) 



This expression contains contributions of one-pair exci- 
tations and those of plasmon excitations. The frequency 
integral in Eq. 1221) can be numerically performed to 
obtain S pp (q). The effective interaction in the RPA, 



(q) is then defined as 



- - 



EE/(p)/(p') 

P,P CTjO-' 



{l-/(p + q)}{l-/(p'-q)> 

£p — e p+q T £p' — e p'— q 



v eff w- 

(23) 



The integral over p and p' in the prefactor has been ana 



lc pr 

lytically performed by Kimball p3 • With the interaction 
v eff A (q) thus calculated we can give a new expression for 
the RPA correlation energy as 



E 



RPA 



1 / 1 



2 \n J 

q.p.p' t7.tr 



E/(p)/(p') 



x Kq) { 1 -/(p+q)Hi-/(p--q)}, eT/(q) . 



£p+q i £p' e p'— q 



(24) 
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FIG. 3: Goldstone energy diagrams with an infinite series of 
particle-partricle ladder interactions and their exchange coun- 
terparts. 



Note that the exchange term is missing in the RPA. The 
ratio of the effective interaction v^jTf A (q) to the bare 
Coulomb interaction v(q) is quadratic for small q, in- 
creases monotonically with q and tends to unity for large 

q- 

In Sec. Ill C and D we employ the expression for 
5g (q) given by Ea. (|22[) as a starting point to allow 
for short-range correlation and exchange. 



Interpolation between long- and short-range 
correlations 



As has been mentioned in Sec. Ill B, the effective 



interaction evaluated in the RPA, 



RPA 
J eff 



(q) is reduced 



in magnitude from the bare Coulomb interaction w(q) 
on the small side of wavenumbers under the influence of 
dynamical screening in the RPA. In fact, the effective 
interaction v e ff(q) is also reduced in magnitude from 
v (q) on the large side of wavenumbers under the strong 
influence of short-range correlation at metallic densities; 
it is not so drastic for small q but essential for large q. 

Consider those perturbation terms which consist only 
of an infinite series of particle-particle ladder interactions 
[5lll52l.l53l . l54ll55j | via the bare Coulomb interaction. Fig. 
3 gives a diagrammatic representation of those pertur- 
bation terms. The correlation energy contribution from 
those diagrams is written as follows: 



E. 



ladd(d) 



= £ £/(p)/(p') 

q,p,p' cr, a' 

.{l-/(p + q)}{l-/(p'-q)} 



2 \n 

x v(q) 



where the particle-particle ladder interaction 7(p,p';q) 
is the solution of the following integral equation: 



£p ^p+q ' ^-p' 



c p -q 



^(p,p';q) = w(q) 
k 



-Tv(a- .^-/(p + kmf-Ap'-k)} 



k 

x /(p,p';k). 



£ p+k 



(27) 



The spin-antiparallel and spin-parallel pair correlation 



functions, .g^(r) and g^(r) derived from E, 

T^laddfex) . . . n 

hie are given respectively as 



ladd(d) 



and 



9i 



r L« = (^) ££/(p)/(p') 



p p' 



1 _{l-/(p + q)}{l-/(p'-q)} 



1 + 



Xl(p,p';q)e' q - r 



(28) 



^) 2 ££/(p)/(p')^|(i- e l(p - p ' 



p p 



1 {l-/(p + q)}{l-/( p '-q)} 



e P e p+q 



c p -q 



I(P, P'; q)e 4q r (l - e l (p-p'+2q) r) _ ^ 



^(p,p';q), 

(25) 



jjjladd(ex) 



--(- 



x u(q) 



£ £<W(p)/(p') 

q,p,p' <j,(t' 

{l-/(p + q)}{l-/(p'-q)} 



€p ^p+q ' £p f ^p ; — q 

x i"(p,p';-p + p'-q). 



(26) 



Eq. (|29|l includes the Hartree-Fock contribution. Both of 

9ladd( r ) and 9H F ( r ) + &9iadd( r )> so far as their behaviors 
at short distances are concerned, are physically sound 
and reasonable throughout the entire region of metallic 
densities in contrast with the RPA in which both of the 
pair correlation functions become negative. 

Consider the general properties of J(p,p';q). It is 
reduced to the bare Coulomb interaction w(q) in the 
limit of small q. On the other hand, it tends approxi- 
mately to w(q)[gjj d(j (0)] 1 / 2 in the limit of large q. The 
ratio I(p, p'; q)/u(q) starts from unity at q = and de- 
creases monotonically with increasing q. Thus, /(p, p'; q) 
is reduced in magnitude from u(q) on the large side of 
wavenumbers under the strong influence of short-range 
correlation at metallic densities. The inclusion of an infi- 
nite series of particle-particle ladder interactions has the 
effect to reduce the contribution from large wavenum- 
ber components of one of the two Coulomb interactions 
entering the second-order direct perturbation term. 

A very good analytic expression for gj^ dd (r) at short 
distances was given by one of the authors (H.Y.) [5l| . 
The magnitude of gj^ dd (0) is remarkably reduced from 
unity as the electron density is lowered to metallic lev- 
els. The inclusion of particle-particle ladder interactions 
much reduces the ratio u e //(q)/w(q) from unity on the 
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large side of wavenumbers, compared with the case of the 
RPA. 

It should be indicated that among all many-body per- 
turbation terms only particle-particle ladder interacting 
parts (via the bare Coulomb interaction) can make a con- 
tribution to the predominant asymptotic form of S'Ti(q) 
of order q for large q. The predominant asymptotic 
form of 5j*(q), on the other hand, is of order q~ 6 for 
large q as a consequence of the cancellation between a 
direct and exchange pair of particle-particle ladder inter- 
actions. 

Consider the iterative solution of the integral equation 
for J(p.p';q). The n-th iterative term involves n + 1 
different wavenumber components of the Coulomb inter- 
actions. Each wavenumber can vary in magnitude from 
zero to infinity under the condition that the total sum 
of ri + 1 different wavenumbers amounts to q. In this 
sense the mixing of different wavenumber components in 
-f(p>P';q) is most complete among all types of interac- 
tions. It is this property of the particle-particle ladder 
interactions that is responsible for the description of local 
density fluctuations. The particle-particle ladder approx- 
imation is in marked contrast with the RPA appropriate 
for the description of long-range correlation associated 
with small wavenumber components of the Coulomb in- 
teraction. The RPA amounts to an infinite summation 
of the most divergent terms of each order in the per- 
turbation expansion and repeats infinite times the same 
wavenumber component q of the Coulomb interaction in- 
stead of mixing different wave number components. 

Let us approximate I(p, p';q) as u(q)(l — C(q)) by 
averaging over p and p' within the Fermi sphere. The lo- 
cal field factor from short-range Coulomb repulsion C(q) 
thus defined is quadratic for small q, increases monotoni- 
cally with q and tends approximately to a constant value 

of i - WiLm 1 ' 2 - 

We make a smooth interpolation between long-range 
correlation in the RPA and short-range correlation in the 
particle-partricle ladder approximation by replacing ev- 
ery bare Coulomb interaction w(q) appearing in Eq.J23) 
ty«(q)(l-C(q))as 



duj- 



1 



at r s — 2 and 0.326 at r s = 4, in the particle-particle 
ladder approximation |5l|. 

The exchange counterpart of particle-particle ladder 
interactions shown in Fig. 3 will be considered in Sec. Ill 
D, combined with other exchange counterparts associated 
with small and intermediate wavenumbers. 



D. Allowance for the exchange correction 

In this subsection we allow for the exchange correc- 
tion, starting with the interpolated effective interaction 
we have obtained above. First of all, we would like to re- 
mind the reader of the detailed analysis of the dielectric 
formulation we have made in terms of the exchange and 
direct parts of the proper polarization function, 7r a (q, uj) 
and 7Tf,(q, uj) in Sec. II. B. So far we have treated the 
Lindhard function 7To(q, lu) that is the non-interacting 
contribution of 7r a (q, uj), the main contribution of the di- 
rect part 7Tf, (q, uj) involving an infinite series of particle- 
particle ladder interactions and the corresponding im- 
proper polarization function, leaving the interacting ex- 
change part 7r a (q, uj) — 7To(q, uj) out of account. In Sec. 
II. B we have stressed that the exchange part 7r a (q, uj) 
has an influence on S'^(q) and 5^(q) indirectly through 
the improper polarization function common to the two. 
This means that the two functions S'^(q) and .^(q) 
are interconnected through ir a (q, uj) in the improper po- 
larization function and have to be determined in a self- 
consistent way. For convenience we term the requirement 
of this self-consistency through the improper polarization 
function the indirect, or feedback effect of exchange. 

In order to allow for the indirect effect of exchange 
as well as the direct one we employ an iterative 
method. First, we start with an exchange interac- 
tion of v e j f RPA+iadd(d) _|_ p / which is obtained 
by replacing q in Ve ^R pA + ladd { d ) (q) by the variable 

P + p' Q appropriate for the exchange term. Sec- 
ondly, we define the local field factor from exchange G(q) 
as follows: 



1 



x tan 



v(q){l - C(q)} 
Hq)(l-C(q))Imh(q^)] 
l + U (q)(l-C(q))Re[7ro(q,c)- 1 > 



1 

n 



E /(p)/(p')<w 

/(p + q)}{l-/(p'-q)} 



It is evident from the interpolation that this expression is 
reduced to Sf? PA (q) for small q and tends to 5c" (q) 
for large q. Following the same procedure as in the RPA, 
we have evaluated the effective interaction appropriate 
for the interpolation. The ratio of the interpolated effec- 

RPA+laddld) , \ , / \ ■ j , r 

(q) to v(q) is quadratic for 



— e 

2 



p+q 



Ep'-q 



v e ff(~p + p' - q) 



E 

p,er,p',( 



/(p)/(p')<W 



{l-/(p + q)}{l-/(p'-q)} 



«(q)G(q). 



(31) 



tive interaction v 



eff 



increases mono- 

Tl (c\X\l/2 



small q in a similar way as in the RPA 
tonically with q and tends to a constant of [^'^(O)] 
which is remarkably smaller than unity at metallic den- 
sities. The value of [g/i^O)] 1 / 2 is estimated to be 0.547 



This is the first value of S^ x (q) to be used in the iter- 
ative process. To make the numerical calculation of the 
left-hand side much easier we have employed the same 
transformation of the variables that Geldart and Taylor 
|56| made in their study of the electron liquid. 
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Thirdly, we replace every factor of u(q)((l — C(q)) ap- 
pearing in Eq.® by u(q)(l - G(q) - G(q)/2) to allow 
for the feedback effect of exchange; a factor of 1/2 is to be 
attached because the local field correction from exchange 
occurs between electrons with the same spin orientations 
alone. The resulting expression is given as 



gRPA+ladd(d)+ex /a _ _g ffF (q) 

n r°° . i 

H 

irn 



duj- 



x tan 



v( q )(l-C(q)-G(q)/2) 
_! »(q)(l - C(q) - G(q)/2)Im[7To(q, «)] 



l + «(q)(l-C(q) 



G(q)/2)Re[7ro(q,w)]' 
(32) 



Fourthly, we evaluate a new effective interaction 
V ™(q) from the direct part of s* PA+ladd{d)+ex {q), 

(q) - B? (q) as 



i.e., 5, 



i?.PA+iadd(d)+e2: 



"n U 



£ /(p)/(p') 



{i-/(p + q)}{i-/(p , -q)l zvw/ , 



(33) 



The interaction (q) thus defined includes the feed- 



back effect of exchange. From v^j j w (q) we construct the 
exchange interaction v^ w (— p + p' q) to calculate a 
new value of the local field factor from exchange G(q). 
We have repeated the iterative process above until we 
reach the convergent result of the local field factor G(q) 
and the effective interaction v e f / (q) ■ 

It should be noted that the above self-consistent de- 
termination of the exchange correction corresponds to 
seeking an approximate solution to the integral equa- 
tion for the particle-hole irreducible interaction within 
the present form of the correlation energy. 

The ratio of the self-consistently determined effective 
interaction u e //(q) to the bare Coulomb interaction w(q) 
is shown for various values of r s in Fig. 4. A remarkable 
change in the ratio w e //(q)/w(q) from its interpolated 
value in Sec. Ill C can be seen over a wide range of 
intermediate wavenumbers. The feedback effect of ex- 
change makes the ratio w e //(q)/u(q) have a broad peak 
around q/pF = 1-4, though its limiting value for large q 
remains unchanged. 

The difference of the ratio v e //(q)/w(q) from unity can 
be interpreted to be a measure of how third- and higher- 
order perturbation terms effectively reduce in magnitude 
one of the two Coulomb interactions entering the second- 
order direct perturbation term. The feedback effect of 
exchange enhances the effective interaction over a wide 
range of intermediate wavenumbers, particularly around 
q/pF = 1-4 almost independently of r s throughout the 
entire region of metallic densities. The origin of the broad 




FIG. 4: The ratio of the effective interaction to the bare 
Coulomb interaction, v e ff (q)A?(q) as a function of q/pF for 
various values of r s . 



peak around q/pF = 1-4 can be traced back to the in- 
teracting exchange part 7r a (q, lu) — 7i"o(q, cj). This can 
easily be seen from the starting expression for the func- 
tion S n (q) in Sec. II B. 



s Ti (q) 



■> ' lm \~r~ : :[ 1 — r~~r ) ~ ^(q,^) 

1 h f°° 

= o<S c (q) - o — / (Mm(7r a (q,tj) - 7r Q (q,w)) 

2 27m J 

= i(S c (q)-5 c (q)). (34) 



The broad peak in the ratio w e //(q)/w(q) comes from 

the peak of S c (q) or S* c (q) located at the almost same 
wavenumber. 



E. Coupling-constant-averaged spin-parallel and 
spin-antiparallel pair correlation functions 

An exact expression for the exchange and correla- 
tion energy functional in DFT is given in terms of 
the cou plin g-constant-averaged pair correlation function 
<?(r,r') |57l l58l 15^. l60| whose behaviors are analogous to 
those of the usual spin-averaged pair correlation function 
ff(r.r')- 



^ 1 f , , ,n(r)n(r') r , .. 

E xc = -jdrdr> ^_^ {ff(r,r')-l}. 



(35) 



This is the inhomogeneous version of the formula for the 
exchange and correlation energy of the uniform electron 
liquid. Similarly, the integral over the coupling constant 



implies the addition of the correlational increase in the ki- 
netic energy to the correlational lowering in the potential 
energy between electrons. The function g(r, r') is defined 
such that the presence of a coupling-constant dependent 
fictitious one-electron potential introduced in a start- 
ing reference Hamiltonian maintains the electron density 
n(r) of the real interacting system while the Coulomb 
interaction among electrons is adiabatically switched on 
as a perturbation. This artifice reflects that correlation 
among electrons in inhomogeneous systems necessarily 
causes a change in the electron density n(r) in contrast 
with the uniform electron liquid. 

The magnitude of the pair correlation function g(r, r') 
at short separations depends on whether one-electron or- 
bitals concerned are extended or localized. Probably, the 
lowering of g(r, r') from 1/2 at short separations will be 
greater for nearly free electrons than for atomic electrons 
tightly bound to the nucleus. This property will also 
apply to the function g(r, r'). 

In this subsection we shall first discuss what effect the 
correlation energy functional generally has on the elec- 
tron density n(r) of non-uniform electron systems with 
the help of the expression for E xc above. The sum of 
the Hartree energy functional and the above E xc gives 
the total electron-electron interaction energy functional 
in the ground state energy functional in DFT. The gen- 
eral properties of g(r, r') are quite analogous to those of 
the usual spin-averaged pair correlation function g(r, r'), 
though the quantity g(r,r') — 1 is somewhat reduced in 
magnitude under the influence of the correlational in- 
crease in the kinetic energy. 

In the Hartree-Fock approximation, g(r,r) = 1/2 and 
1/2 < g(r, r') < 1. Generally, correlation reduces the 
value of g(r,r') on the short-distance side of |r — r'| 
and reversely enhances the value of g(r,r r ) on the long- 
distance side such that the total electron charge involved 
is conserved; at zero separation < g(r, r) < 1/2 and 
the function g(r, r') increases with |r — r'|, exceeds unity 
and forms a broad peak around the average inter-electron 
separation and thereafter asymptotically approaches to 
unity. The function g(r, r') can be rewritten as g(f, R) 
with r = r r' and R = (r + r') /2 and its value at zero 
separation <?(0, R) should vary rather slowly from place 
to place in solids. 

Repelling of other electrons from the immediate neigh- 
borhood of an electron in question, when it approaches 
each nucleus in solids, makes it "feel" the attractive po- 
tential from the nucleus significantly stronger because of 
less screening than in the Hartree-Fock approximaton in 
DFT. Hence the electron density n(r) is enhanced in the 
immediate neighborhood of the nuclei. Generally, cor- 
relation causes an enhancement in the electron density 
n(r) at the positions of individual nuclei in solids. 

Next we shall mention the general properties of 
the coupling-constant-averaged spin-parallel and spin- 
antiparallel pair correlation functions for the uniform 
electron liquid, g''(r) and g'^{r) and their Fourier trans- 
forms iS^(q) and <S^(q) before we show their numerical 
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results calculated from the present theory. The coupling- 
constant-averaged pair correlation functions g aa (r) are 
defined as 

i r 2 

g™ (r) = - / d(e 2 )g™ (r). (36) 
e Jo 

The integral over the coupling constant generally re- 
duces both the shift of g^ (r) from unity and the shift 
of g^(r) from its Hatree-Fock value over all distances 
since the correlational increase in the kinetic energy 
is involved. The reduction is conspicuous particularly 
at short distances. The cusp condition on g 00 (r) 

mm mm mum is modmed as 




The difference between g 00 (r) and g aa (r) in the cusp 
condition is the presence of the coupling constant integral 
operators in the case of g aa (r) . The cusp conditions on 
g^(r) and g^(r) are fulfilled in the particle-particle lad- 
der approximation and in the Hartree-Fock plus paricle- 
particle ladder approximation, respectively. We have al- 
ready given the expression for g^(r) and g'^(r) in these 
approximations in Eas. l|28|) and (|29|l . respectively. These 
cusp conditions hold also for inhomogeneous electron sys- 
tems if one averages the value of g aa (r, r') on the sphere 
centered at R = (r + r')/2 with a small radius |r — r'| 

[HI- 

In Fig. 5 the value of g^ (r) calculated from the present 
theory is drawn for various values of r s . We have ascer- 
tained that the value of g^(r) is significantly enhanced 
at short distances, compared with the value of <?^(r) cal- 
culated in the particle-particle ladder approximation. As 
is obvious from Fig. 6, the present g*' (r) fulfills the Pauli 
principle. We have also observed that the shift of g^ (r) 
from its Hartree-Fock value is reduced at short and inter- 
mediate distances, compared with the case of g^(r) cal- 
culated in the same approximation. Perdew and Wang 
m| have also calculated the spin-averaged g(r) of the 
electron liquid by interpolating between the RPA and the 
particle-particle ladder approximation, but they have left 
the interacting exchange part 7r Q (q, ui) — 7To(q, ui) out of 
account. 

Fig. 7 shows the functions ST (q) and SJJ-(q) in the 
present theory. For comparison the same functions in 
the RPA are also drawn in the figure. These functions 
are both reduced in magnitude from the corresponding 
functions without bars for the same reason as in the case 
of g aa '(r) L The difference between ^ T (q) and 5^(q) is 

given by S c (q). For lack of the second- and higher-order 
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FIG. 5: The coupling-constant-averaged spin-antiparallel 
pair correlatrion function for the electron liquid, <?^(r) calcu- 
lated from the present theory as a function of pft for various 
values of r s . 
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FIG. 7: (a) and (b) The coupling-constant-averaged spin- 
parallel and spin-antiparallel correlational contributions of 
the static structure factor, 31 (q) and 5c (q) as functions 
of q/pF for r 3 = 2.0 and 4.0; the dotted line in each figure 
represents the value of the two functions in the RPA. 



FIG. 6: The coupling-constant-averaged spin-parallel pair 
correlatrion function for the electron liquid, <? TT (r) calculated 
from the present theory as a function of pft for various values 
of r s . 



exchange corrections the RPA fails to discriminate be- 
tween the two functions. As is obvious from Fig. 7, the 
inclusion of exchange corrections shifts the two functions 
in opposite directions for small q. This splitting is evi- 
dent from the very definition of S'^(q) and St-'-(q) that 
we have given in Sec. II B. The interacting exchange 
part of proper polarization function 7r a (q, u>) — 7To(q, uj) 
is responsible for this splitting of 5*J T (q) and S'l^q) for 
small q. In this connection we mention that the accurate 
description of S*(q) in the RPA for small q does not imply 
that it also can give the accurate description of S(q) for 



small q. 

The function <SP(q) becomes positive for q > l.bpF 
while the function iS^q) stays negative for any q. This 
is because the exchange contribution overcomes in mag- 
nitude the direct contribution for intermediate and large 
wavenumbers. The following asymptotic forms of S'J^q) 
and i5J*(q) for large q are an alternative representation 
of the cusp conditions [55| on <7^(r) and (r). 

^W-|(?)(?)V(P) + .., 09) 

ft)' « 

It is straightfoward to derive the corresponding asymp- 
totic forms of §1^ (q) and Si' (q) for large q. 
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FIG. 8: The function 5J* (q) for various values of r s 
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FIG. 9: The function Sj^q) for various values of ; 



In Figs. 8 and 9 5j^(q) and «§,F(q) are drawn for var- 
ious values of r s , respectively. The correlation energy of 
the electron liquid can be evaluated by integrating 5^(q) 
and iSJ^(q) over q, or equivalently by the area of the two 
curves drawn in Figs. 8 and 9 relative to the horizontal 
axis. The value of the correlation energy calculated from 
the present theory agrees with the most reliable values 
within an accuracy of 0.5 mRy. per electron through- 
out the entire region of metallic densities (see Table I). 
However, one must not forget that the accuracy of the 
correlation energy is not necessarily the best criterion for 
the validity of a theory. For example, the self-consistent 
GW approximation gives a very accurate evaluation of 
the correlation energy of the electron liquid, but it makes 
the value of <7(r) negative at short distances 68]. 

The present theory gives an accurate description of 
both S'j^q) and 5jf(q), or equivalently their Fourier 
transforms. Furthermore, it gives with high accuracy 
spin-parallel and spin-antiparallel components of the cor- 
relation energy e c over the entire region of metallic den- 
sities. Fig. 10 shows the ratios of e c **/e c and ej' /e c as 
functions of r s . In the high density region, the two ratios 
are analytically given as 

e c 2 2 0.06221n(r s ) " ' ' 

e c TT _ 1 1 c 2x 
e c ~ 2 2 0.0622 ln(r s ) 
e c = 0.0622 ln(r s ) + c RPA + c 2x + ■ ■ ■ Ry., 
CRPA = -0.1422, c 2x = 0.04836. (41) 

where crpa is the constant from the RPA and c 2x is the 
value of the second-order exchange energy [6|| ■ It is im- 
portant to recognize that e]} and eJ T occupy nearly 70% 
and 30% of the total correlation energy throughout al- 
most the entire region of metallic densities, respectively; 
these ratios probably will apply to all valence electrons 



°- 8 | i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 
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r s 

FIG. 10: The ratios of t]} /e c and eV 1 /e c as functions of r s . 



participating in cohesion of solids. The component e\} of 
e c comes from the deepened Fermi hole at short and in- 
termediate distances from its Hartree-Fock value, as can 
be seen from g*'(r). The component e^, on the other 
hand, comes from the developed Coulomb hole. The par- 
allel and antiparallel components equivalently contribute 
to the leading logarithmic term in the high density limit. 
This is because iSj-'-(q) and Sl'(q} become equivalent in 
the limit r s — » 0. This tendency can be observed from 
the two functions drawn in Figs. 8 and 9 for r s = 0.1. 
In the language of the pair correlation functions, the 
leading logarithmic term comes from such distances as 
1 <C PFf *C r s where g^{r) — 1 and <?^(r) — g]j F (r) 
both behave like —r s /(pFr) 2 ; the screening length is pro- 
portional to r s /pf for small r s |55| . 
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r a 


RPA 


STLS 


C-A 


VWN 


EPX 


Present 


1.0 


-157.6 


-123 


-119.6 


-120.6 


-119 


-120.7 


2.0 


-123.6 


-91.4 


-90.2 


-89.6 


-89.1 


-89.9 


3.0 


-105.5 


-74.8 




-73.8 


-73.7 


-73.7 


4.0 


-93.6 


-64.0 




-63.6 


-63.6 


-63.3 


5.0 


-85.0 


-56.3 


-56.3 


-56.3 


-56.3 


-55.8 


6.0 


-78.2 


-50.5 




-50.7 


-50.7 


-50.1 



TABLE I: Values of the correlation energy calculated from a 
variety of approximations and methods for the electron liquid, 
RPA, STLS 0, C-A, VWN, EPX and the present interpo- 
lation method; C-A and VWN stand for the original values 
Ceperley and Alder obtained from the Green's function Monte 
Carlo method (GFMC) and the values from an interpolation 
formula using their GFMC data, respectively; EPX stands for 
the effective potential expansion method proposed by Takada 




The approximation method of replacing g(r, r') in the 
exact formula for E xc in Eq. (|35fl by its analogue of 
the uniform electron liquid has been pro posed by a num- 
ber of authors [7J El El El El III This method 
has the merit of giving the nonlocal effect of E xc in the 
sense that the resulting v xc (r) implicitly depends on the 
electron density n(r') at the neighboring position r' of 
r. As for the density parameter to be used in the func- 
tion g(r — r') of the electron liquid, we think that the 
uniform electron density in the interstitial region out- 
side muffin-tin spheres for metals or the averaged den- 
sity of the system as a whole is rather better than the 
local density. In this case it will be appropriate to ap- 
proximate the functional derivative Sg(v' — r")/<5n(r) as 
S(t - (r' + r") /2)dg{r' - r")/dn where dg{r' - r")/dn 
should be evaluated at the same density parameter as 
in g(r - r'). 

F. Some criticisms on the interpolation 



and made a smooth extrapolation to smaller wavenum- 
bers assuming an appropriate form quadratic in the limit 

There remains some arbitrariness in the choice of the 
local field factor C(q) around intermediate wavenumbers. 
Those proper polarization functions of the direct type 
77b (q,cj) which involve second- and higher-order particle- 
particle ladder interactions make the predominant contri- 
bution to S'^(q) and S^(q) for intermediate and large 
wavenumbers. It should, however, be pointed out that 
those proper polarization functions of the same direct 
type which involve particle-hole and hole-hole ladder in- 
teractions can make an appreciable contribution for in- 
termediate wavenumbers (see (3h) and (3j) in Fig. 1). 
By its strict definition, the local field factor from short- 
range Coulomb repulsion C(q) should be determined 
from all the proper polarization functions of the direct 
type 7Tfe(q,w). 

Since the first attempt by Hubbard [lil Il5| , a number 
of interpolation methods have been proposed with the lo- 
cal field correction to the RPA for the study of the elec- 
tron liquid in the region of metallic densities. The present 
theory is the most sophisticated among these attempts. 
Here we have managed to circumvent the treatment of 
muti-pair excitations |22l |23| by evaluating the coupling- 
constant-averaged static structure factor S(q) instead of 
the static structure factor 5*(q). If one starts with S(q) in 
his interpolation, one cannot avoid treating multi-pair ex- 
citations and thereafter has to perform the integral over 
the coupling constant to enumerate correctly the many- 
body perturbation terms involving particle-particle lad- 
der interactions. Instead, we have started with the exact 
ground state energy contribution from an infinite series of 
particle-particle ladder diagrams in order to interpolate 
between the RPA and the ladder approximation in the 
form of the coupling-constant-averaged static structure 
factor S(q). 



Finally, we offer some criticisms of the present interpo- 
lation method for its possible development in the future. 
We have first interpolated between long-range correlation 
in the RPA and short-range correlation in the particle- 
particle ladder approximation in Sec. Ill C and then 
allowed self-consistently for the exchange correction to- 
gether with its feedback effect on the direct interaction 
by iteration in Sec. Ill D. Once the local field factor 
from short-range Coulomb repulsion C(q) introduced in 
Sec. Ill C is given, there is no problem in the present 
method of evaluating the exchange correction; the local 
field factor from exchange G(q) is uniquely determined. 

The present theory depends on the accuracy of the 
starting interpolation, i.e., the local field factor from 
short-range Coulomb repulsion C (q) . One of the authors 
(H.Y.) [51j gave an analytic solution to the ladder inter- 
action I(p, p'; q) averaged over p and p' within the Fermi 
sphere. We have employed this analytic solution for the 
determination of the factor C(q) primarily for q > pp 



IV. ORBITAL-DEPENDENT CORRELATION 
ENERGY FUNCTIONAL 

A. Orbital-dependent exchange energy functional 

According to the optimized-effective-potential method 
(OEP) |8( in DFT, the ground state energy functional of 
many-electron systems can be written as 

N r / fi 2 \ 

E a = £ /drtf (r H__ V 2 W(r) 

2 — 1 ^ ' 

+ J dvv ext (r)n{v) + \jdv J dr'^^1 

2jjJ J |r-r'| 

+ E c [{^}}. (42) 
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The first term is the kinetic energy functional of the 
reference noninteracting system with the same electron 
density n(r) as the real system; this term is usually de- 
noted by T s . The second is the electron- nucleus inter- 
action energy functional. The third and the fourth are 
the Hartree and the exchange energy functionals, respec- 
tively. The exchange energy functional ^[{^j}] is rep- 
resented as an explicit functional of Kohn-Sham orbitals 
and hence £7 x [{y?j}] is an implicit functional of the elec- 
tron density n(r); i and j have the same spin orienta- 
tions. The sum of the third and fourth terms amounts to 
the electron-electron interaction energy functional in the 
Hartree-Fock approximation. We assume that the cor- 
relation energy functional is represented as an explicit 
functional of Kohn-Sham orbitals. It is important to rec- 
ognize that the correlation energy functional E c [{(fi}] de- 
fined in DFT amounts to the net correlational lowering 
in the electron-electron interaction energy partially can- 
celled by the correlational increase in the kinetic energy, 
T c ; E c = T c + V^ e where V e c e is the correlational lowering 
in the electron-electron interaction energy. 

The Hohenberg-Kohn variational principle 
5E g [n(r)]/6n(r) — can be expressed in terms of 
a self-consistent set of Kohn-Sham equations describing 
the reference noninteracting system as 



(43) 



- — V 2 + v s {r) } ipi{r) = e m (r), 
2m ' 



v s (r) = v ext (r) + v H (r) + v x (r) + v c (r) , (44) 



v H (r) = J dr' 



where v s (r) denotes the total single electron potential; 
v e xt(r) is the Coulomb potential from the nuclei and 
vh{y) is the Hartree potential; v x (r) and v c (r) denote 
the exchange and the correlation potentials defined by 
v x (r) = SE x [{ipi}]/5n(r) and v c (r) = 5E c [{<pi}]/8n(r), 
respectively. 

In this subsection we mention only the exchange energy 
functional E x [{ipi}] and its functional derivative v x (r) 
since it is a starting point fundamental to the introduc- 
tion of an orbital-dependent correlation energy functional 
in DFT. An excellent review on the theoretical treatment 
of the orbital-dependent exchange energy functional in 
DFT and its numerical results for atoms, molecules and 
solids has recently been given by Grabo, Kreibich, Kurth 
and Gross [8j. 

If one employs the orbital-dependent exchange energy 
functional, one has to evaluate v x (r) as follows: 



^(r) = 5E x [{(fi}]/Sn(r) 

- ST f ^^r j ^E x [{ Vi }]5 Vi {v') 
~ 2fJ I <V*(r') Mr") 



where 6v s (r')/5n(r) is the inverse of the static density 
response function of the reference noninteracting system, 
X _1 (r, r'). The functional derivative Sipi(r')/5v s (r") and 
the response function x( r > r ') can both be expressed in 
terms of Kohn-Sham orbitals and Kohn-Sham energies. 
After a few manipulations one can obtain the integral 
equation which determines the exchange potential v x (r). 

OCC. r, 

Y, j dx 1 {ip\{Y')G sl {Y' ,v) Vl (v) + C .c}v x {y') 



Sv a (r") 
Sn(r) 

(45) 



V /rfr'l^fr'l — ^Mfi ( r > 



r)<^( r ) 



G si (r,r')=£ 



1j( r )y]( r ') 



+ c.c > , (46) 



(47) 



It is straightfoward to generalize this integral equation 
to include the correlation potential v c (r) as well if an 
orbital-dependent correlation energy functional £J C [{^}] 
is available. 

An alternative method can lead to the same integral 
equation for the exchange-correlation potential v xc (r)(— 
v x (r) + v c (r)). Consider the ground state energy func- 
tional Eg in Eq. I|42|l under the condition that (pi (r) 's are 
the solutions of Eg. (|43fl with an unknown one-electron 
potential v s (r) . Then, E g can be considered to be a func- 
tional of v s (r) and there arises a variational problem of 
how to optimize v s (r) in the minimization of E g , i.e., 
SEg[v s (r)]/Sv s (r) = 0. The potential v s (r) thus calcu- 
lated is named the optimized-effective-potential (OEP). 
From this optimization one can also derive the same inte- 
gral equation that we have just mentioned above. Hence 
it is termed the OEP integral equation. In fact, the 
optimization SE g [v s (r)]/Sv s (r) — is equivalent to the 
Hohenberg-Kohn variational principle 5E g [n(r)]/6n(r) = 
|77l l78| . We summarize the conclusions obtained from 
the so-called exchange-only OEP method in the follow- 
ing. 

(1) The exchange-only OEP method gives the exact 
exchange potential v x (r) that completely cancels the self- 
interaction terms involved in the Hartree potential. 

(2) The derivative discontinuity in ^[{t^}], or the dis- 
continuity of v x (r) can be correctly described. 

(3) The Hartree-Fock ground state energy in the frame- 
work of DFT is by only a slight amount higher than the 
usual Hartree-Fock ground state energy. This is because 
in DFT the additional condition of a local form is im- 
posed on the exchange potential in the minimization of 
the same ground state energy functional. 

(4) Instead, the electron density n(r) calculated from 
DFT might be by only a slight amount better than the 
one calculated from the usual Hartree-Fock theory owing 
to the variational principle. 

However, it is important to recognize that the 
exchange-only OEP integral equation is identical to the 
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first iterative form of the exchange-only Sham-Schluter 
equation (see Eq. l|58l) '). The Sham-Schluter equation im- 
plies that the standard many-body theory based on the 
one-electron Green's function formalism and DFT both 
lead to the same exact electron density n(r) of many- 
electron systems. This requirement of course should be 
satisfied only when both of the two many-body theories 
are exact. If applied to the Hartree-Fock stage of ap- 
proximation, this equation amounts to forcibly requiring 
that the Hartree-Fock approximation in the framework of 
DFT lead to the same electron density n(r) as the usual 
Hartree-Fock theory; it cannot be fulfilled. From the first 
iterative form of the exchange-only Sham-Schluter equa- 
tion it is evident that the difference in the electron den- 
sity n(r) between the two methods will be only a slight 
amount just like the difference in t he g round state energy. 

(5) Li, Krieger and Iafrate [7!j have applied the 
exchange-only spin-unrestricted OEP method to atoms 
(Z=3-56(Li-Ba), 79 and 80 (Au and Hg)) to make 
a comparison with the spin-unrestricted Hartree-Fock 
method (SUHF). Their calculations are very accurate; 
the virial theorem and the related identity (V* e + 
J drn(r)r ■ gradual - ) = 0) for the exchange energy func- 
tional E x (= V* e ) and the exchange potential v x (r) due 
to Levy and Perdew [8(| are satisfied to four or five sig- 
nificant figures. 

The ground state energy calculated from the exchange- 
only OEP method is at most by 40ppm higher than 
that from the Hartree-Fock method and the differnce 
decreases with increasing Z. They report that the elec- 
tron density n(r) calculated from the exchange-only spin- 
unrestricted OEP method is close to the SUHF value. 
Interestingly, they observe a few reversals of the order 
of the state yielding the highest occupied eigenvalue be- 
tween the ns and (n— l)d states in the transition series of 
atoms. The Kohn-Sham exchange potential v x {v) always 
raises the energy levels of the inner shell states compared 
with those of the SUHF; the same trend can also be seen 
from a comparison between the Hartree-Fock energy lev- 
els of occupied states of the electron liquid and the corre- 
sponding Kohn-Sham energy levels. The reversals occur 
only when the ns and (n— l)d energy levels are very close 
in the SUHF. Along the spirit of the Slater approximation 
to the exchange potential, Krieger et al (KLI) j8ll . 183. l83| 
have proposed a simple approximate solution to the OEP 
integral equation in which the complexity of the original 
OEP integral equation is significantly reduced but many 
of its essential properties are kept unchanged. 

(6) Buijse and Baerends [84| write that the Hartree- 
Fock theory significantly distorts the electron density 
distribution n(r) of atoms to the extent that often the 
largest absolute errors in the total Hartree-Fock energy 
occur for the electron-nucleus interaction energy and the 
kinetic energy rather than for the electron-electron inter- 
action energy and that the Slater approximation for the 



exchange potential gives rather better results for the elec- 
tron density n(r) of atoms. However, the largest relative 
error, i.e., the largest ratio of each error to the corre- 
sponding Hartree-Fock value probably will occur for the 
electron-electron interaction energy. 

(7) Hartree-Fock calculations make higher 4s or- 
bitals occupy prior to lower 3d orbitals in 3d transi- 
tion atoms. This reversal of the occupied states has 
been supposed to be ascribed to the variational principle 
5E^ F [{(fii}]/S(pi = 0. The reversal has been interpreted 
to be a consequence of the strategy that the total energy 
will not be raised due to the electron-electron interaction 
energy in E^ F [{tfi}]. On the other hand, no reversal 
of occupied states should occur if one accurately solves 
the exchange-only OEP integral equation in DFT. This is 
because Kohn-Sham equations by definition are intended 
for the reference non-interacting system. This problem 
may be closely related to the central field approximation 
conventionally applied to atoms. However, no systematic 
method has yet been developed to solve the OEP inte- 
gral equation appropriately for a nonspherically symmet- 
ric system. Strictly, the degeneracy of 3d or 4d orbitals 
should be lifted for lack of spherical symmetry and ns 
and (n—l)d orbitals should be mixed as a result. 

(8) As for the energy gap of semiconductors, the 
exchange-only OEP method gives a good result [83, ISg . 

(9) Strictly, the Hartree-Fock theory and its density- 
functional version by the x-only OEP method generally 
give different results of the ground state energy E g , the 
electron density n(r) and the ionization potential ( or the 
chemical potential ); for the uniform electron liquid the 
two give the same result of these quantities. 



B. Correlation-induced reconstruction of 
Kohn-Sham orbitals 

In this subsection we propose an explicitly orbital- 
dependent correlation energy functional beyond the 
Hartree-Fock approximation for the study of the elec- 
tronic structure in atoms, molecules and solids in DFT. 
It is to be employed in combination with the orbital- 
dependent exchange energy functional. 

In Sec. Ill we have defined a new expression for the 
correlation energy of the electron liquid. It is a modifi- 
cation of second-order perturbation terms, in which one 
of the two Coulomb interactions in each term is replaced 
by an effective interaction containing information about 
long-, intermediate-, and short-range correlations. Sim- 
ilarly, we prop ose here a second-order perturbation like 
expression 87] which is constructed from Kohn-Sham or- 
bitals, Kohn-Sham energies, and the effective interaction 
borrowed from the electron liquid. 
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1 occ. unocc. s 



{LPiip 3 \v{ri 2 )\ip a 'Pb){^a'Pb\veff{r 1 2)\'Pi'p J ) (iptcpj \v(ri 2 ) \iPaVb) (Patpb\v e ff (^12) I WPi) 



i,j a,b 



e i + e j — e a — e b £j + e j — e a — e 6 J 

(48) 

{(pi(Pj\v(r 12 )\(pa<Pb) = / dridr293*( r i)^( r 2)w(ri2)^Q(ri)(pb(r 2 ), (49) 



(ViVjlwe//^)!^^} = / dr 1 dr2ip*(r 1 )(p* j (r 2 )v e ff(r 1 2)(pa{ri)(pb(r2), (50) 



where i,j and a, 6 represent occupied and unoccupied 
Kohn-Sham states, respectively and the states i,j in the 
second-order exchange like term have the same spin ori- 
entations; v(r\2) is the Coulomb interaction for r\2 — 
n — r 2| The effective interaction w e //(?'i2) in Ea. lf^ 
is the real-space Fourier transform of the interaction 
v e ff(q) we have evaluated for the electron liquid in 
Sec. III. How to choose in practical calculations the den- 
sity parameter for the effective interaction we borrow 
from the electron liquid will be discussed in detail in Sec. 
IV C. 

Here we note that a fully orbital-dependent correlation 
energy functional, or equivalently an orbital-dependent 
effective interaction in the present scheme can in principle 
be obtained if one follows the same procedure as in the 
electron liquid with Kohn-Sham orbitals and Kohn-Sham 
energies in place of planewaves and free electron energies. 

From a comparison with a general expression for 
the exchange-correlation energy functional using the 
coupling-constant-averaged pair correlation function 
g(r, r') (see Eq.JSSJl) it may be naturally expected that 
the present correlation energy functional has also the ef- 
fect to reduce the Hartree-Fock electron-electron interac- 
tion energy in its contributions from short distances. Let 
us characterize the present orbital-dependent correlation 
energy functional. 

(1) The exact many-electron wavefunction fulfills the 
requirement of antisymmetry. In accordance with this 
antisymmetry, the electron-electron interaction energy 
functional in DFT should be constructed from direct and 
exchange pairs of terms. The present orbital-dependent 
correlation energy functional as well as a direct and ex- 
change pair of electron-electron interaction energy func- 
tionals in the Hartree-Fock approximation reflects the re- 
quirement of antisymmetry. 

(2) In the limit of uniform density the present the- 
ory is reduced to the local density approximation (LDA) 
as it should be. The self-interaction terms involved in 
the Hartree energy are exactly cancelled by the orbital- 
dependent exchange energy functional. Furthermore, 
higher-order self-interaction terms are, from the begin- 
ning, eliminated from the present -E c [{<p.;}] since it con- 
sists of a direct and exchange pair of terms. Therefore the 
present theory is completely free from self-interaction er- 
rors. This is indispensable for the treatment of exchange 



and correlation between electrons with localized orbitals. 
The LDA estimates the location of the occupied 3d-bands 
in copper and zinc too highly because of self-interaction 
errors. The overestimation as large as about 18% in the 
neck cross section of the Fermi surface of copper in the 
[111] direction, which is considered to be ascribed to the 
too high 3d-bands in the LDA, will be for the most part 
improved by removing self- interaction errors |8Sj . 

(3) One of the remarkable merits of the present the- 
ory is that the physical meaning of correlation is clarified 
by the present representation of the functional £? c [{<^j}]. 
According to the virial theorem, the correlation energy of 
atoms, molecules and solids under no constraint of any 
external force amounts to a half of the lowering in the 
total potential energy since it is by half cancelled by the 
correlational increase in the kinetic energy. In this sense 
the partitioning of the ground state energy functional 
in DFT is helpful in understanding the meaning of the 
correlation energy since it is an exact many-body theory 
assuming the self-consistent one-electron formalism. The 
total potential energy consists of three terms. The first 
is the electron-electron interaction energy and the second 
the electron-nucleus interaction energy. The third is the 
nucleus-nucleus interaction energy. Usually, the first and 
second potential energies are both lowered at the cost of 
increasing the corresponding two kinetic energies, T c and 
T s . Here we assume the Born-Oppenheimer approxima- 
tion. 

It is important to notice that correlation usually 
changes the magnitude of the nucleus- nucleus interac- 
tion energy from its Hartree-Fock value. The Hartree- 
Fock approximation regularly underestimates the equi- 
librium length between the nuclei in molecules. The LDA 
that borrows knowledge of the exchange-correlation en- 
ergy from the uniform electron liquid has a tendency to 
underestimate the Wigner-Seitz radius of metals or the 
lattice constant of solids. Probably, these shortcomings 
will be closely related to the fact that both approxima- 
tions overestimate the screening of the nuclei at short 
distances. 

By definition, the functional -^[{t^}] in DFT repre- 
sents the net lowering of the electron-electron interaction 
energy partially cancelled by the correlational increase 
in the kinetic energy T c . From a comparison with the 
Hartree-Fock total energy functional it is evident that 
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the direct term in Ea. (|48(l has the effect to reduce the 
Hartree energy and its exchange counterpart has the ef- 
fect to reduce in magnitude the exchange energy in a 
well-balanced way through the same effective interac- 
tion v e ff(r) that contains accurate information about 
long-, intermediate-, and short-range correlations beyond 
second-order perturbation theory. 

Furthermore, it is important to recognize that the pres- 
ence of -E c [{</?i}] in the total energy functional not only 
has the effect to reduce the net electron-electron interac- 
tion energy in its contributions from short distances but 
also causes a change in the electron density n(r) such 
that the total energy is further stabilized by a lower- 
ing in the only negative energy, i.e., the electron-nucleus 
interaction energy though it inevitably increases the ki- 
netic energy of the reference non- interacting system T s . 
The partitioning of the ground state energy functional in 
DFT has the merit of connecting the net lowering in the 
electron-electron interaction energy (in its contributions 
from short distances) with the further lowering in the 
electron-nucleus interaction energy through Kohn-Sham 
equations that minimize the total energy functional. 

(4) The present theory is expected to give an accurate 
evaluation of the correlation energy even of the most lo- 
calized electrons such as atomic electrons in which the 
LDA makes the most serious errors because of its failure 
to cancel self-interaction terms and its overestimation of 
the correlation energy. Then we may expect that the 
present B c [{^}] is valid for tightly binding electrons as 
well as for nearly free electrons. These two different types 
of electrons are both participating in the cohesion of tran- 
sition metals. An accurate evaluation of the ground state 
energy of atoms, molecules and solids may be generally 
expected from the present theory if the density parameter 
for the effective interaction borrowed from the electron 
liquid is suitably chosen. The optimization of the density 
parameter will be discussed in detail in Sec. IV C. 

(5) One of the greatest merits of DFT lies in the 
fact that one can describe all correlation effects on one- 
electron states with the local correlation potential v c (r) 
in the framework of self-consistent one-electron theory, 
without being troubled with the treatment of the non- 
local, energy-dependent, and complex self-energy opera- 
tor involving many-body background excitations. This 
merit is acquired only when accurate information about 
the potential n c (r) is available. The total single elec- 
tron potential v s (r) in Kohn-Sham equations determines 
the electron density n(r) of the system. Its accuracy 
depends entirely on the exchange-correlation potential 
v xc (r). The exchange potential v x (r) can be exactly 
treated by the exchange-only OEP method. Consider 
what effect the correlation potential v c (r) has on the elec- 
tron density n(r) of the system. 

The correlated motion of valence electrons in solids, 
though it inevitably increases the kinetic energy, more 
lowers the electron-electron interaction energy in its 
contributions from short distances. The simultaneous 
change in the electron density n(r), i.e., the shrinkage of 



the electron distribution around the nuclei arising from 
less screening of the nuclei due to short-range correlation, 
further lowers the only negative electron-nucleus interac- 
tion energy at the cost of increasing the kinetic energy of 
the reference non- interacting system T s . Then the low- 
est possible ground state energy is realized such that the 
nucleus-nucleus interaction energy is lowered by some- 
what increasing the lattice constant as a counteraction 
of the shrinkage of the electron distribution around the 
nuclei and the increase in T s is to some extent relaxed 
accordingly. 

Thus, correlation in real solids is necessarily accompa- 
nied by a change in the electron density n(r) in contrast 
with the case of the uniform electron liquid. It may be 
generally expected that correlation induces inhomogene- 
ity in the electron density distribution n(r) of the system 
against the Hartree-Fock approximation in DFT. 

In the neighborhood of individual nuclei of the sys- 
tem, the total single electron potential v s (r) in Kohn- 
Sham equations can be written in terms of the coupling- 
constant-averaged pair correlation function <?(r,r') as 
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The screening of individual nuclei at short-, intermediate- 
, and long-distances can be properly described by the 
second and third terms; for isolated systems like atoms 
these two terms give the correct asymptotic form v s (r) = 
— e 2 /|r| for large |r| since the following sum rule holds: 



dr'n{r'){g(r,r') - 1} = -1. 



(52) 



The last term in Eq. (|51|l probably will not give a signifi- 
cant r-dependence for any |r| and plays a background-like 
role in v s (r) since the integrand including the functional 
derivative Sg(r' ,r")/Sn(r) is integrated over r' and r". 
Consider the physical role of <?(r,r') in the sum of the 
second and third terms when the distance |r — r' | is de- 
creased. Since the inequality that < ,g(r,r') < 1/2 
holds for |r — r'| less than the average inter-electron sep- 
aration, it is evident that the Hartree or Hartree-Fock 
screening of individual nuclei is significantly reduced at 
short distances due to the presence of <?(r, r'). In other 
words, short-range correlation yields less screening of the 
attractive nuclear potential at short distances. It is this 
less screening of the nuclear potential that shrinks the 
electron distribution in the vicinity of the nuclei to en- 
hance the electron density n(r) at the positions of indi- 
vidual nuclei. 

If r and r' are both in the interstitial region outside 
muffin-tin spheres of metals, the function g(r,r') is very 
well approximated by the corresponding function of the 
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uniform electron liquid we have calculated in Sec. Ill E 
and hence causes no change in the electron density. On 
the other hand, if r and r' are both inside a muffin-tin 
sphere, the probability of g(r,r') is concerned with not 
only valence electrons but also core electrons. The mag- 
nitude of g(r, r') at short separations within the muffin- 
tin sphere is less than 1/2 due to short-range correla- 
tion, though it may be somewhat greater than in the 
interstitial region. Short-range correlation in this essen- 
tially nonuniform density region reduces the Hartree or 
Hartree-Fock screening of the nuclei and hence causes the 
shrinkage of the electron density distribution around the 
nuclei to lower the electron-nucleus interaction energy 
more than its Hartree-Fock value. 

(6) Most importantly, the present orbital-dependent 
correlation energy functional i£ c [{y>j}] hybridizes 
Hartree-Fock Kohn-Sham orbitals in order to produce 
sclf-consistently reconstructed Kohn-Sham orbitals 
through its functional derivative 5E c [{ipi}]/5n(r), 
i.e., the correlation potential v c (r). Naturally, this 
correlation-induced reconstruction of Kohn-Sham or- 
bitals gives a further stabilization of the ground state 
energy as well as a redistribution of the electron density 
n(r) of the system. The present correlation potential 
v c (r) is expected to produce less screening of the 
attractive nuclear potential at short distances to give 
a shrinkage of n(r) around individual nuclei, thereby 
leading to a lowering in the electron-nucleus interaction 
energy. 

The Hartree-Fock approximation, in its original ver- 
sion or its DFT version, overestimates the effect of the 
electron-electron interaction on the screening of the nu- 
clei at short distances as well as the magnitude of the 
electron-electron interaction energy. Variation of the ex- 
pectation value of the total Hamiltonian within a sin- 
gle Slater determinant overestimates the influence of 
the interaction on occupied orbitals in the form of self- 
consistent Coulomb and (nonlocal or local) exchange po- 
tentials. The Hartree-Fock orbitals are then spread too 
far away from the nuclei and hence the kinetic and the 
electron-nucleus interaction energies in this approxima- 
tion are both underestimated in magnitude. 

From these peculiarities it is evident that this approxi- 
mation tends to underestimate nucleus-nucleus separa- 
tions in order to get a sufficient lowering in the only 
negative electron-nucleus interaction energy for the sta- 
bilization of molecules. There are some cases such as 
F2, in which the Hartree-Fock approximation is unable 
to stabilize molecules and the redistribution of the elec- 
tron density due to correlation is indispensable to the 
formation of molecules. 

The LDA is exact in the limit of uniform density but, 
when applied to nonuniform electronic systems, it cannot 
describe properly such exchange-correlation effects as are 
peculiar to nonuniform systems. The most important of 
these effects is correlation-induced less screening of the 
nuclei at short distances. The LDA tends to underesti- 
mate the lattice constant of solids for lack of this effect 



and furthermore tends to overestimate the magnitude of 
the cohesive energy because of a rather too uniform spa- 
tial variation of the exchange-correlation potential v xc (r). 
Notice that the LDA considerably smooths out the spa- 
tial variation of v xc (r) since it borrows the knowledge of 
fi xc (n(r)) from the uniform electron liquid; fi xc {n(r)) is 
the exchange-correlation contribution of the chemical po- 
tential of the uniform electron liquid appropriate for the 
electron density n(r). 

(7) On the basis of the virial theorem and his own 
partitioning of the total energy, Ruedenberg has gained 
a profound insight into the binding nature of molecules 
[8!j,|9(|. He has emphasized that the binding of molecules 
or the cohesion of solids (including metals) is essentially 
stabilized by shrinking the electron density distribution 
in the immediate vicinity of the nuclei to enhance the 
electron density n(r) at the positions of indivisual nuclei, 
compared with the case of isolated atoms. 

According to Moruzzi, Janak and Williams [7J, on the 
other hand, the electron density in simple metals is re- 
versely reduced in the immediate vicinity of the nuclei. 
This is because the LDA fails to give the proper descrip- 
tion of the correlation-induced less screening of individ- 
ual nuclei at short distances. In the LDA the Hartree 
screening of the nuclear potential is simply shifted by 
an amount of fi xc (n(r)). The density dependence of 
fjt X c(n(v)) appropriate for the electron liquid is abso- 
lutely insufficient for describing the correlation-induced 
less screening of nuclei. It is beyond the scope of the 
LDA since the electron density n(r) in the vicinity of the 
nuclei in metals is far from being slowly varying. 

The present theory is promising for the realization of 
the enhancement of the electron density in the vicinity 
of the nuclei as a consequence of less screening of the 
nuclear potential. This is because the present £J C [{^}] 
involves all unoccupied excited states that are needed to 
redistribute the electron density n(r) in this region; the 
resulting v c (r) is expected to have much larger negative 
values than the LDA correlation potential in this region. 

Kohn-Sham one-electron states have to be recon- 
structed while atoms condense into solids. We think that 
the reconstruction of Kohn-Sham onc-clcctron states as- 
sociated with the cohesion of solids can also be well de- 
scribed by the present correlation energy functional pro- 
vided the density parameter in the effective potential is 
appropriately chosen. 

In connection with the considerations above we men- 
tion a notewo rthy fact. According to Boyd and his 
coworkers [tnl l92t 193$ , Hund's multiplicity rule can be 
ascribed to a gain in the electron-nucleus potential en- 
ergy caused by correlation-induced less screening of the 
nuclei in contrast with the traditionally accepted inter- 
pretation due to the exchange energy. On the contrary, 
the electron-electron interaction energy is grater in a high 
spin state for all neutral atoms and molecules calculated 
with requisite accuracy. They have arrived at this con- 
clusion from accurate configuration interaction calcula- 
tions of low-lying excitation energies of light atoms and 
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molecules, in which the virial theorem is of course satis- 
fied. 

(8) There is a strong possibility that the electronic 
structure of the so-called strongly correlated electron sys- 
tems may be interpreted in terms of the correlation- 
induced reconstruction of Kohn-Sham orbitals. The 
energy-band structure of solids calculated in the LDA 
is approximately equivalent to the band-structure that 
is calculated in the Hartree approximation and shifted 
to the low energy side by an amount of fj, xc of the 
uniform electron liquid with the averaged electron den- 
sity of solids. In marked contrast to the LDA, the 
present orbital-dependent correlation energy functional 
hybridizes Hartree-Fock Kohn-Sham one-electron states 
with different orbital characters, particularly in the 
neighborhood of the Fermi level or the energy gap, in 
order to produce self-consistently reconstructed Kohn- 
Sham one-electron states. The occurrence of such recon- 
struction of one-electron states is peculiar to non- uniform 
real systems, quite absent in the uniform electron liquid 
which consists of electrons with the same character. 

Consider the energy denominator in the present 
E c [{<fi}}. Evidently, the functional l? c [{y>j}] can gain its 
magnitude if the density of states is sharply enhanced in 
the immediate neighborhood of the Fermi level under the 
influence of self-consistently reconstructed Kohn-Sham 
energies. Such an enhancement in the density of states 
should be accompanied by a considerable change in the 
electron density n(r) particularly for strongly anisotropic 
systems. 

Our scenario of heavy fermions is as follows: As a 
consequence of strong hybridization between conduction 
electrons and f electrons, the occupied and unoccupied 
states will probably be reconstructed such that they en- 
hance the density of states in the immediate neighbor- 
hood of the Fermi level as sharply as possible in order to 
gain the magnitude of the functional 2? c [{^}], which is 
balanced against a simultaneous increase in the kinetic 
energy of the reference non-interacting system. 

For the study of the electronic structure of heavy 
fermion systems it is of course indispensable to adopt 
the fully relativistic one-electron Dirac equations |94| in 
the framework of DFT. 

(9) Almbladh and von Barth [95| have studied the re- 
lations between Kohn-Sham eigenvalues and exact exci- 
tation energies. From exact asymptotic results for the 
spin density far away from a finite system, they have 
shown that the uppermost Kohn-Sham eigenvalue equals 
the ionization potential, not only for extended systems 
like solids, but also for finite systems like atoms. 

From available accurate wavefunctions for a number of 
light atoms (the He isoelectronic series, Li and Be), Alm- 
bladh and Pedroza 96] have constructed E xc , v xc (r) and 
v c (r) to evaluate the Kohn-Sham eigenvalues. They have 
then found that the uppermost (or highest occupied) 
Kohn-Sham eigenvalue reproduces the exact ionization 
potential and that the eigenvalues for deeper shells lie 
above the corresponding exact excitation energies, while 



the Hartree-Fock eigenvalues lie below. 

The correlation potential v c (r) they have constructed 
for He, Li + and Bc ++ takes large negative values in the 
immediate vicinity of the nucleus reflecting that the exact 
electron density is slightly higher than the Hartree-Fock 
value in this region and on the other hand becomes posi- 
tive at large r from the nucleus, forms a broad peak and 
tends to zero with increasing r. This is quite different 
from the potential v c (r) in the LDA since it is always 
negative and increases monotonically with increasing r. 

Generally, the Koopmans theorem gives a fairly good 
evaluation of the ionization potential for atoms and 
molecules, though Hartree-Fock eigenvalues take no ac- 
count of both correlation and relaxation ^jfy- This is 
because a great amount of cancellation occurs between 
the extra correlation energy of the ./V-body ground state 
relative to the (N — l)-body excited state and the relax- 
ation energy (which is also a kind of correlation energy) 
of the (N — l)-body excited state produced by removing 
an electron from the Hartree-Fock N-body ground state. 

This cancellation is not complete, however, and the 
Koopmans theorem usually overestimates the ionization 
potential, or equivalently the Hartree-Fock theory tends 
to evaluate the highest occupied eigenvalue (negative) 
too low. This implies that the relaxation energy is usu- 
ally somewhat larger than the extra correlation energy in 
magnitude; for much localized orbitals like the He atom it 
is so large that the Koopmans theorem cannot be applied. 
Thus, the correlation correction to the ionization poten- 
tial is usually negative for molecules; the Hartree-Fock 
Kohn-Sham equations with the optimized local exchange 
potential will lead to the same conclusion. The presence 
of nuclei in atoms and molecules is responsible for the 
occurrence of relaxation. If an electron is removed, the 
remaining N — 1 electrons have to readjust themselves to 
the attractive nuclei. 

It is very reasonable that the correlation potential v c (r) 
for He, Li + and Be ++ takes positive values at large dis- 
tances from the nucleus if one considers that the expec- 
tation value of i> c (r) with respect to the highest occupied 
state has to give a correction with the positive sign. 

On the other hand, the solutions of Hartree-Fock or 
Hartree-Fock Kohn-Sham equations for the uniform elec- 
tron liquid are plane waves and hence there occur no 
relaxation of orbitals. Then, only the extra correlation 
energy is missing in this case. The correlation contribu- 
tion to the chemical potential, /i c for the uniform electron 
liquid is thus negative. @@@@@@@@ @ 

Averill and Painter [9^, and Levy and Perdew |80| 
have given an expression for the virial theorem in terms 
of density-functional theory. Their expression for the 
virial theorem at the equilibrium positions of the nuclei 
is splitted into two parts, i.e., the reference one-electron 
(non-interacting) and the correlational parts. 

2T C + V e c e + J dm(r)r • gradu c (r) = 0. (53) 

This is the correlational part of the virial theorem that 
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relates the correlational increase in the kinetic energy 
T c , the correlational lowering in the electron-electron 
repulsion energy V e c e , and the correlation potential v c (r); 
T c is usually incorporated into the exchange-correlation 
energy functional E XC (E XC = T c + V^ e ). This equation is 
helpful to discuss the general shape of v c (r). 

Let us first consider the quantity 2T C + V^ e in Eci. (|53[l . 
It vanishes to second order in e 2 for atoms and molecules. 
It also vanishes for the electron liquid in the high den- 
sity region where the correlation energy is given by the 
leading logarithmic and the subsequent constant terms. 
Consider the contribution to 2T C + V£ e from very large 
momentum transfer interactions in the electron liquid 
at metallic densities, or very short-range parts of the 
Coulomb interaction. It also vanishes since the relation 
2T C + V e c e = holds for this case 0. The reason for 
this cancellation is that the particle-particle ladder in- 
teraction I(p, p';q), or multiple scattering between two 
particle states with opposite spins makes a contribution 
to both T c and V e c e in a way analogous to second-order 
perturbation theory. This type of cancellation between 
2T C and V e c e , associated with short-range correlation, ap- 
plies also to atoms, molecules and solids. 

From these considerations, 2T C + V e c e is a rather small 
but negative quantity. It is negative for the electron liq- 
uid at metallic densities. Second-order perturbation the- 
ory for atoms and molecules, by which 2T C + V^ e vanishes, 
tends to overestimate the magnitude of the correlation 
energy and hence the total higher-order contribution to 
2T C + V° e is expected to be negative. We may then con- 
clude from Eq.J^Sl) that the integral J dvn{r)v ■ gradi> c (r) 
is a small positive quantity for atoms, molecules and 
solids (including metals). 

Correlation changes the electron density n(r) in atoms, 
molecules and solids. The correlation potential v c (r) is 
responsible for this change in n(r). Short-range corre- 
lation rectifies the overestimated screening of the nuclei 
in the Hartree- or Hartree-Fock approximation and en- 
hances the electron density n(r) in the immediate vicin- 
ity of the nuclei as a result. Correspondingly, the po- 
tential v c (r) takes large negative values in this region. 
As the distance from the nucleus r increases, the value 
of u c (r) increases with a steep positive slope. In order 
to conserve total electronic charge, the potential v c (r) 
for an atom has to become positive for larger r from 
the nucleus, forms a broad peak and tends to zero with 
increasing r. The correlation potential v c (r) that Alm- 
bladh and Pedroza have constructed from a few accurate 
many-electron wavefunctions actually shows such a non- 
monotonic behavior. 

Similarly, the correlation potential v c (r) for metals will 
probably take large negative values in the immediate 
vicinity of each nucleus and steeply increase with the dis- 
tance from the nucleus. It is expected to exceed ndK) 
(instead of zero) within the muffin-tin sphere as the dis- 
tance from the nucleus r is increased, form a broad peak 
and tend to /x c (r*) as r approaches the muffin-tin sphere 
radius; /x c (r*) is the correlation contribution of the chem- 



ical potential of the uniform electron liquid with the den- 
sity parameter r* appropriate for the interstitial region 
of metals. It should be noted that the potential v c (r) in 
the immediate vicinity of nuclei probably will take larger 
negative values for metals than for isolated atoms, re- 
flecting that the electron density in this region should be 
enhanced under the influence of cohesion. 

The tendency of the LDA to overestimate the cohesive 
energy of metals in magnitude and to underestimate the 
lattice constant of metals provably will be ascribed to 
the failure of its correlation energy functional to give the 
above nonmonotonic behavior of the correlation poten- 
tial v c (r) as well as to the failure of its exchange energy 
functional to cancel spurious self- interaction terms. 

(10) The correlation energy defined in DFT, E® FT is 
in magnitude slightly larger than the usual correlation 
energy, E c defined by the differnce between the usual 
Hartree-Fock and the exact ground state energies (see 
IV A. (3)). According to the virial theorem, E^ F1 - = 
-{T S -T? F )-T C = i(V es + V£ + V* e )-±(V es + V*r™ F 
where V es is the total classical electrostatic energy that 
consists of the nucleus-electron, the Hartree, and the 
nucleus-nucleus interaction energies; V* e and V e c e are the 
exchange and correlation contributions to the electron- 
electron interaction energy, respectively. The inequalities 
that T s > T HF > T" F hold since two different virial the- 
orems with the same formal expression are satisfied for 
the usual Hartree-Fock approximation and its density- 
functional version (see Ea l54|) . The quantity T HF stands 
for the kinetic energy of the usual Hartree-Fock theory. 

We give an explanation of the correlation energy de- 
fined in DFT, E° FT . The first term, -(T s - T S HF ) is 
the negative of an increase in the kinetic energy of the 
reference non-interacting system, which is primarily as- 
sociated with a correlation-induced enhancement in the 
electron density in the immediate vicinity of nuclei, giv- 
ing a correlational lowering in the nucleus-electron and 
nucleus-nucleus interaction energies. The second term, 
— T c , on the other hand, is the negative of an increase 
in the kinetic energy arising from the correlated motion 
of many electrons, which is primarily associated with a 
correlational lowering in the electron-electron interaction 
energy. 



C. Density parameter appropriate for the effective 
interaction 

In Sec. Ill B we have mentioned that a fully orbital- 
dependent correlation energy functional, or an orbital- 
dependent effective interaction can in principle be ob- 
tained if one makes a similar interpolation between long- 
and short-range correlations involving exchange effects 
with Kohn-Sham orbitals and Kohn-Sham energies in 
place of planewaves and free-electron energies for the 
electron liquid. At present, such a sophisticated inter- 
polation is exceedingly difficult to perform. Instead, we 
borrow the knowledge of the effective interaction from 
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the electron liquid. In this case there arises a problem 
of how to choose the density parameter for the effective 
interaction. We must then resolve this problem to make 
the present theory applicable to practical calculations. 
For this purpose we adopt the following method. 

Moruzzi, Janak and Williams Q were the first who 
successfully applied DFT to energy-band calculations of 
metals with the atomic number less than 50 using the 
LDA. One of the most important conclusions they have 
obtained from their systematic study is that the com- 
pressibility of the electron liquid calculated for the uni- 
form electron density in the interstitial region outside 
muffin-tin spheres is generally a good approximation to 
the compressibility of real metals. This general trend 
can be observed over a wide variety of metals covering 
simple, transition and noble metals. The value of the 
inverse compressibility 1/k for transition metals is gen- 
erally large owing to the presence of d electrons, while the 
value of 1/k for alkali metals is small. It is very small 
particularly for sodium, potassium, and rubidium with 
the lowest metallic densities. 

We may therefore conclude that the electron liquid 
model is valid for all metals in the sense that the uni- 
form electron density in the interstitial region makes a 
predominant contribution to the inverse compressibility 
1/k of metals. This implies that the metal cannot survive 
if its uniform electron density in the interstitial region ex- 
ceeds the lowest possible critical density of the electron 
liquid where the compressibility k becomes divergent and 
the system can no longer be stable thermodynamically. 
This critical density corresponds to about r s = 5.25. It 
has been generally observed that the uniform electron 
density in the interstitial region of metals is somewhat 
enhanced in comparison with the averaged density of va- 
lence electrons as a whole. Let us describe the uniform 
electron density in the interstitial region by the effective 
density parameter r s *. Then, r s * < r s for almost all 
metals. The values of r s * for sodium, potassium, and ru- 
bidium are about 3.66, 4.36, and 4.67, respectively, while 
those of r s are 3.96, 4.87, and 5.18, respectively. For 
cesium with the lowest observable metallic density, the 
usual density parameter r s = 5.57 is beyond 5.25, but its 
effective density parameter r s * is considered to be less 
than 5.25. 

In fact, the LDA has borrowed the knowledge of the 
exchange-correlation energy from the electron liquid be- 
yond its applicability. One can practically calculate its 
ground state energy beyond the critical density resorting 
to a variety of approximation methods, but this does not 
implies that the system stays stable. According to many- 
body theoretic techniques, the exact theory should main- 
tain the internal consistency between self-energy and ver- 
tex parts and fulfill the Pauli principle. Since the vertex 
part is closely related to the compressibility k of the sys- 
tem, this internal consistency can no longer be main- 
tained when it is divergent. In this sense the electron 
liquid is no longer stable for r s > 5.25. 

In accordance with the situation above, we borrow the 



knowledge of the effective interaction from the electron 
liquid for r s < 5.25. A practical choice of the density 
parameter for the effective interaction probably will be 
the parameter r* appropriate for the uniform electron 
density in the interstitial region of metals calculated from 
the LDA. The best choice of the density parameter is 
the optimization that minimizes the calculated ground 
state energy in the present theory as a function of the 
density parameter; r* will be a good approximation to 
the optimized density parameter. 

The present correlation energy functional borrows the 
effective interaction from the electron liquid. There- 
fore we have to treat separately the functional derivative 
5v e ff(r')/Sn(r) when we solve the OEP integral equa- 
tion in order to evaluate the exchange-correlation poten- 
tial v xc (r)(= v x (r) + v c (r)). An approximation appro- 
priate for this purpose will be that Sv e f /(r')/<5n(r) = 
S(r — r')dv e ff(r)/dn. The derivative dv e f f(r)/dn should 
be evaluated at the same density parameter as in the 
effective interaction v e ff(r). 

So far we have discussed how to choose the density 
parameter for the effective interaction v e ff(r) when it is 
intended for use in metals. Originally, the knowledge of 
the correlation energy of the electron liquid should be 
applied to the study of the electronic structure of met- 
als. The LDA, if applied to atoms, underestimates the 
exchange energy in magnitude by 14 ~ 8% and reversely 
over estim ates the correlation energy by a factor of 2 ~ 2.5 
[99l llOOj . This is because the spurious self-interaction 
terms are not cancelled and because the correlation en- 
ergy of the uniform electron liquid increases its magni- 
tude like lnr s Ry. per electron at high densities. 

In reality, the correlation energy of atoms is in mag- 
nitude comparable with that of conduction electrons in 
metals, i.e., about 75 ~ 66 mRy. per electron. However, 
the ratio of the correlation energy to the exchange energy 
is an order of magnitude larger in conduction electrons 
than in atomic electrons; in atoms the accurate evalua- 
tion of the exchange energy is indispensable before the 
correlation energy is allowed for. 

The present theory gives the exact evaluation of the 
exchange energy and moreover does not overestimate the 
correlation energy of atoms at all since it borrows the 
knowledge of the effective interaction from the electron 
liquid, not of the density dependence of the correlation 
energy. It should, however, be noted that the effec- 
tive interaction v e ff(r) of the electron liquid, if applied 
to atoms, probably will lead to an underestimation at 
long distances. This is because the screening in atoms 
is not perfect in contrast to the electron liquid. If one 
defines the effective interaction properly for atoms in a 
similar way as in the electron liquid, the ratio of the 
Fourier transformed effective interaction v e ff(q) to the 
bare Coulomb should tend to a finite value of order 1/eo, 
not to zero in the limit q — > where eo is the static 
dielectric constant of atoms. Fortunately, the small- 
wavenumber components of i; e //(q) are not expected to 
make an important contribution to the correlation energy 
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of atoms since the predominant part of the correlation 
energy comes from virtual transitions to unoccupied lo- 
calized orbitals with discrete spectra. Similarly, the ratio 
of the effective interaction properly defined for insulators 
and semiconductors should have a finite value in the limit 
q — > owing to the presence of the energy gap. 

On the other hand, no qualitative difference in the 
intermediate- and short-wavenumber components of the 
effective interaction may be expected between the elec- 
tron liquid and real atoms, insulators, and semiconduc- 
tors. 

To conclude, the best choice of the density parameter 
for the effective interaction in the present theory is the 
optimization that minimizes the calculated ground state 
energy as a function of the density parameter, irrespec- 
tive of the difference between of metals, insulators and 
atoms. It is noteworthy that in the limit of uniform den- 
sity Kohn-Sham orbitals are reduced to plane waves and 
the resulting ground state energy in the present theory 
takes its minimum when the density parameter for the 
effective interaction amounts approximately to r s = 4.2. 

There is another method of choosing the density pa- 
rameter for v e f f (r) . It utilizes the reference one-electron 
part of the virial theorem in terms of DFT |80l l9Sj . 

2T S + V es + V* - J drn(r)r • gradt> c (r) = 0, (54) 

where V es (V es — V ne + Vj^ + V nn ) is the total classi- 
cal electrostatic energy functional that consists of the 
nucleus-electron interaction energy V ne , the Hartree en- 
ergy Vg^f and the nucleus- nucleus interaction energy V nn ; 
Vgg stands for the exchange energy functional (V^ e = E x ). 
The density parameter for v e f / (r) can in principle be de- 
termined from the requirement that Eq. l|54|l be satisfied. 
The minimization of the exact ground state energy ex- 
pression and the virial theorem should be fulfilled simul- 
taneously at the exact equilibrium inter-nucleus separa- 
tion. 

It should be noticed that the usual expression for the 
virial theorem, 2T + V — can be obtained from Eqs. 
(UHl and g2J; T = T s + T c and V stands for the total 
potential energy of the system. 

In the electron liquid with the neutralizing uniform 
positive background correlation shifts the equilibrium 
density parameter from r s — 4.83 to 4.2. Correspond- 
ingly, the ground state energy is lowered from —0.095 
to about —0.154 Ry. per electron. From the virial 
theorem with and without correlation it follows that 
2(T - T HF ) + - V£ F ) = where T and V ee are 
the exact kinetic and electron-electron interaction ener- 
gies of the electron liquid, respectively; the superscript 
HF is attached to the corresponding Hartree-Fock val- 
ues. 

There is a second-order perturbation like relation [54| 
between large wave number components of T associated 
with the one-electron momentum distribution function 
n(k) and those of V ee associated with the static structure 



factor 5(k) due to short-range correlation, i.e. 2T(k) + 
F ee (k) = 0for |k| >p F . 

In real molecules and solids, on the other hand, the 
Hartree-Fock model generally gives a fairly good approx- 
imation to the equilibrium nucleus-nucleus separation. 
Correlation tends to increase the equilibrium separation 
by only a small amount. The Hartree-Fock potential 
curve of diatomic molecules as a function of the nucleus- 
nucleus separation, if shifted nearly uniformly downward 
by an amount of the correlation energy, leads to the ex- 
act potential curve in the neighborhood of the equilib- 
rium separation. This implies that 2(T — T HF ) + (V ee — 
V e H e F ) + (V ne - V n H /) (or equivalently -(V nn - V n H n F ) ) 
nearly vanishes or is a small positive value where V ne 
and V nn denote the electron-nucleus and nucleus-nucleus 
interaction energies, respectively. 

Besides the net correlational lowering in the electron- 
electron interaction energy partially cancelled by the cor- 
relational increase in the kinetic energy, the correlation- 
induced less screening of the nuclei produces a fur- 
ther lowering in the electron-nucleus potential energy by 
shrinking the electron distribution in the vicinity of the 
nuclei and accordingly the lattice is somewhat expanded 
to reduce the nucleus-nucleus interaction energy, relax- 
ing the simultaneous increase in the kinetic energy of the 
reference non-interacting system T s to some extent. 

Thus, the correlational lowering in the electron- nucleus 
potential energy, which is absent in the uniform elec- 
tron liquid and hence cannot be properly described in 
the LDA, plays a significant role in the accurate evalua- 
tion of the equilibrium nucleus- nucleus separation as well 
as the electron density n(r). 



V. RELATIONSHIP BETWEEN KOHN-SHAM 
EQUATIONS AND DYSON EQUATIONS 

In this section we first evaluate the quasiparticle en- 
ergy dispersion of the electron liquid with high accuracy. 
Next we give a detailed discussion about the relationship 
between Kohn-Sham equations and Dyson equations in 
order to justify the application of Kohn-Sham equations 
to the energy-band theory. 



A. Symmetric expression for the self-energy 

Sham and Kohn were the first who discussed the 
relationship between Kohn-Sham equations and Dyson 
equations. In the system where the electron density is 
slowly varying, they have found that in the neighborhood 
of the Fermi level \x the quasiparticle energy dispersion 
E(p) — [i can be related to the Kohn-Sham energy dis- 
persion e KS (p) — fi as 

E{p)-fi= C -^TW 1 ^ . ( 55 ) 

Jdr2^% p (r)|2' 
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where m*(n(r)) denotes the effective mass of the uniform 
electron liquid appropriate for the local density n(r) and 
ifp (r) is the Kohn-Sham orbital specified by wavenumber 
P 

It is very difficult to give a qualitatively reliable evalu- 
ation of the effective mass to* of the electron liquid over 
the entire region of metallic densities. This is because the 
expansion of the self-energy with respect to the screened 
Coulomb interaction converges very slowly in this region. 
The effective mass ratio to* /to calculated to the first or- 
der in the RPA screened interaction, giving the leading 
ter m in the high density expansion of to* /to for small 
i~ s |lOl| . takes the minimum around r s — 0.7, exceeds 
unity around r s — 2.3 and thereafter slowly increases 
with r s . The ratio to* /to calculated to the second order 
stays less than unity until r s = 6, but it has also the 
minimum around r s = 1 and thereafter increases slowly 
with r s . The very slow convergence of the effective mass 
ratio m* /to reflects that correlation is much complicated 
in the region of metallic densities. 

In the language of many-body perturbation theory, this 
implies that the determination of the quasiparticle en- 
ergy dispersion requires the accurate knowlegde of the 
self-energy E CT (p, e) as a functional of the one-electron 
Green's function G a >(p',e') rather than a few lowest or- 
der expansion of S CT (p, e) in the screened Coulomb inter- 
action, though screening is of course necessary to over- 
come the exchange difficulty. As is proved in a cele- 
brated book by Noziers, any variation in the self-energy 
^cr(Pi e ) occurs only through a variation in each con- 
stituent G CT <(p',e') of the self-energy E CT (p,e). This is 
the reason why the accurate knowledge of S (T (p, e) as a 
functional of G cr /(p / ,e') is essential for the evaluation of 
the quasi-particle energy dispersion. 

Yasuhara and Takada [24| have introduced the very 
expression for the self-energy S CT (p, e) that is symmetric 
with respect to each consitituent G in it. 

a' J y ' J n=l 

x ( P) e;p',e';[G])G CT ,(p',e'), 

OO 

I™' (p, e; p', e'; [G]) = £ i»-'( P) e; p ', e ' ; [ G f]) j 

n=l 

(56) 

where l( n ) crcr (p, e; p', e'; [G]) denotes the contribution 
to the particle-hole irreducible interaction from all 
the possible n-th order skeleton diagrams. Accord- 
ing to this symmetric expression, one can easily check 
that the functional derivative <5E cr (p, e; [G])/5G CT < (p', e') 
is equivalent to the particle-hole irreducible inter- 
action I aa (p, e; p', e'; [G]), as is proved in standard 
textbooks on many-body theory. The interaction 
I aa (p, e; p', e'; [G]) is the value of the full particle-hole ir- 
reducible interaction I aa (p, e; p', e'; q, ui; [G]) in the limit 
of q = lo = 0. The symmetric expression has the strik- 
ing merit of visualizing that a great deal of cancellation 



occurs systematically between higher-order contributions 
due to the Pauli principle. This is indispensable to the 
qualitative and quantitative evaluation of the quasiparti- 
cle energy dispersion as well as the effective mass of the 
electron liquid in the strongly correlated density region. 

A great deal of cancellation can also be observed be- 
tween the wavenumber dependence and the energy de- 
pendence of the self-energy S(p, e; [G]) in the determina- 
tion of the quasiparticle energy dispersion -E'(p). This can 
be traced back to the dynamical nature of screening in 
the electorn liquid. Both of the wavenumber dependence 
and the energy dependence are comparatively large in 
magnitude, but they cancel each other for the most part 
in the determination of E(p) and m*. 

Si nce the ea rly stage of the study of the electron liq- 
uid |l02t Il03j . it has been recognized that the first it- 
erative solution for the quasiparticle energy dispersion, 
i.e., E(p) = e(p) + E CT (p, e(p); [G ]) is a good approxi- 
mation to the solution of the fully self-consistent equa- 
tion, E(p) = e(p) + S ff (p, E(p); [G]). This is because 
there is a great deal of cancellation between the higher- 
order iterative correction and the correction made by 
replacing £<j(p, e; [G ]) with S cr (p, e; [G]); the partially 
self-consistent solution with S CT (p, e; [G ]), the solution 
of e = e(p) + £ CT (p, e; [G ]) becomes worse. Almost no 
difference can be expected between the first iterative and 
the fully self-consistent solutions around the Fermi level 
since the electron liquid is a many-body system which 
consists of electrons with the same character. 

The self- energy as a functional of G°, T, a (p, e; [G ]) can 
be written as [43| 

E CT (p, £ ;[G°]) = ^i |de 2 (e 2 )-t 
xE/p/ d e 7-'(p, £ ;p', e ';[G ])G«,(p', e '). 

(57) 

where the sum over n with the factor l/(2n — 1) in 
Eq. H56|) 15 transformed into the coupling constant in- 
tegral since G° does not depend on the coupling con- 
stant. All diagrams for the particle- hole irreducible in- 
teraction I aa (p, e; p', e'; [G ]) can be classified into three 
different types, imp roper, doubly irreducible, and proper 
and reducible |24| (see Fig. 2). From this classifica- 
tion it is evident that there is a great deal of system- 
atic cancellation among higher order contributions of 
I acF (p, e; p', e'; [G ]) due to the Pauli principle. In the 
practical calculations we have employed an approximate 
form of I aa (p, e; p', e'; [G ]), in which the above cancel- 
lation due to the Pauli principle is allowed for system- 
atically up to higher orders; on the other hand, we have 
put the contribution from J cr_<T (p, e; p', e'; [G ]) as a con- 
stant since it is rather small; we have chosen the constant 
by fitting the value of S(p, e(p); [G ]) at p = pp to the 
exchange-correlation contribution of the chemical poten- 
tial \l X c- 
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FIG. 11: The real part of the self-energy correction calculated 
from the present theory, ReE(p, e(p); [G ]) in units of the 
Fermi energy ef as a function of q/pF for various values of 
r s . For comparison the correction calculated from the GW 
approximation is also drawn. 



Fig. 11 shows the real part of the self-energy correc- 
tion calculated from the present theory in units of the 
Fermi energy for various values of r s . For compari- 
son the correction calculated from the GW approxima- 
tion is also drawn in the figure. A clear difference in 
the wavenumber dependence of ReX^p, e(p); [G ]) can 
be seen between Q < p < pp + q c and pf + q c < P where 
q c is the cutoff wavenumber of plasmon excitations. The 
present self-energy correction Re£ CT (p, e(p); [G ]) main- 
tains a small positive slope over a range ofO < p <PF+q c 
throughout the entire region of metallic densities. The 
slope of ReE CT (p, e(p); [G ]) at the Fermi wavenumber pp 
is closely related to the effective mass to* and is a slowly 
increasing function of r s . The value of the mass ratio 
to* /to estimated from the slope is consistent with the re- 
sult in a previous study of the effective mass 0, 0, EH • 
The mass ratio to* /to is a slowly decreasing function of 
r s and is shifted from unity about 7 ~ 8% even at the 
lowest metallic densities. This means that the wavenum- 
ber dependence of S CT (p, e; [G]) somewhat overcomes in 
magnitude the energy dependence of £ CT (p,e; [G]) in the 
determination of the quasiparticle energy E(p) in the 
neighborhood of the Fermi level. 

For p > pp + q c , on the other hand, the self-energy cor- 
rection ReE (T (p, e(p); [G ]) behaves approximately like 
— 1/p in a quite different way from for < p < pf + q c - 
The critical wavenumber near pp + q c corresponds ap- 
proximately to the energy fj, + huj p i where hw p i is the 
plasmon excitation energy. The remarkable difference in 
the wavenumber dependence of the self-energy correction 
can be ascribed to the fact that screening does work only 
for < p < pf + q c in the electron liquid. 

The GW approximation gives a fairly good evaluation 
of the magnitude of the self-energy correction, but it fails 



to give a qualitatively correct description of its derivative 
for < p < pf + Qc over almost the entire region of metal- 
lic densities. The appearence of a dip around pf + q c in 
both cases is ascribed to the use of Re£ CT (p,e; [G ]) in 
place of Re£ CT (p, E(p); [G]). The dip should be destroyed 
under the influence of a coupling between plasmon exci- 
tations and multipair excitations. 

From these observations we may conclude that the 
self-energy correction to the quasiparticle energy disper- 
sion E(p) is fairly well approximated by its value at 
the Fermi wavenumber pp for a range of wavenumbers 
< p < pf + q c throughout the entire region of metallic 
densities. In the limit of uniform density the Kohn-Sham 
energy dispersion becomes e KS (p) — e(p) + fi xc . The ra- 
tio {ReS(p, e(p); [G°] — //i C }/fe is less than about 10% 
for < p < pp + q c - Therefore we may say that Kohn- 
Sham energies are a rather good approximation to quasi- 
particle energies so far as excitation energies relative to 
the Fermi level are less than about fiw p i . 

The conclusion above will probably apply to real met- 
als. This is because the large cancellation between the 
wavenumber dependence and the energy dependence of 
the self-energy £(p,e;[G]) in the determination of the 
quasiparticle energy dispersion E(p) below the critical 
energy near [i + Hujpi can be traced back to the dynamical 
nature of screening that is associated with the presence 
of the Fermi surface. 

An energy-band calculation of sodium |104| performed 
in the Yasuhara-Takada approximation |24j much more 
advanced than the usual GW approximation has actually 
revealed that the occupied band width of sodium is in- 
creased by an amount less than 10%, compared with the 
conventional LDA result. 

Very recen tly, R . Maezono, M.D. Towler, Y. Lee, and 
R. J. Needs |105| have performed the diffusion Monte 
Carlo (DMC) energy-band calculation of sodium to find 
confidently that its occupied band width is increased to 
the almost same extent as the present result. It seems 
that the exchange difficulty of Hartree-Fock equations 
associated with metals as well as the electron liquid has 
long been believed to be resolved completely by the RPA, 
but in fact it is not the case, i.e., the RPA only removes 
the infinite derivative at pf ■ It is concluded that the 
present many-body theoretic study and the above DMC 
calculation have cooperated to resolve in a truly complete 
sense the exchange difficulty of the Hartree-Fock theory 
at realistic metallic densities. 



B. Relation between £(r, r',e) and v xc (r) 

Dyson equations and Kohn-Sham equations are de- 
rived from two different variational principles concerning 
the same ground state energy E g , SE g /SG(r, r', e) = 
and 5E g /dn(r) — 0, respectively, where E g can be con- 
sidered to be a functional of the one-electron Green's 
functional G(r, r', e) or another functional of the electron 
density n(r). The two different types of equations can in 
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principle give the same exact electron density n(r), chem- 
ical potential and ground state energy of a many-electron 
system. 

Sham and Schluter |l06l Il07j have derived an ex- 
act expression for the exchange-correlation energy func- 
tional E xc by means of many-body theoretic technique 
based on the Green's function formalism and related the 
exchange-correlation potential v xc (r) in DFT to the self- 
energy £(r,r',e; [G]). This Sham-Schluter equation im- 
plies that the electron density n KS (r) evaluated from the 
Kohn-Sham Green's function G KS (r,r',e) as the unper- 
turbated Green's function in the many-body perturba- 
tion expansion is identical to the electron density n(r) 
calculated from the exact one-electron Green's function 
G(r,r',e). 

Grabo, Kreibich, Kurth and Gross |2| have discussed 
the OEP method from the viewpoint of many-body per- 
turbation expansion. They have then pointed out that 
the higher order contribution of the exchange-correlation 
energy functional beyond the Hartree-Fock approxima- 
tion, Ex'c 1 (n ^ 2) depends not only on occupied and 
unoccupied Kohn-Sham orbitals and Kohn-Sham ener- 
gies but also on the exchange-correlation potential v xc (r). 
This implies that the resulting OEP integral equation in- 
volves the functional derivative Sv xc (r)/5n(r') as well as 
the potential v xc (r) and is much more difficult to solve 

as a result. 

■ I 

Casida 108] has started with an approximate expres- 
sion for the ground state energy functional. In the expan- 
sion of the self-energy functional using skeleton diagrams 
he has omitted those terms which arise from the coupling 
constant integral of the constituent renormalized Green's 
functions, or equivalently those terms which vanish if 
the exact one-electron Green's function G is replaced by 
the unperturbed Kohn-Sham Green's function G KS in 
the exact ground state energy functional for the varia- 
tional principle. Using this approximate ground state 
energy functional E [G KS ] + §[G KS ], Casida has gen- 
eralized the exch ange-only OEP method by Sharp and 
Horton [Tol EM 

to include correlation. The resulting 
OEP integral equation has turned out to be identical 
to the first iterative form of the Sham-Schluter equation 
(see Ea. (|58J) '). Casida has thus shown that the exchange- 
correlation potential v xc (r) in DFT can be interpreted 
to be the variationally best optimized local approxima- 
tion to the exchange-correlation part of the self-energy 
£(r,r',e;[G]). 

We would like to stress here that Casida's expression 
for the ground state energy functional Eq[G k s ] + $[G ks ] 
is a sufficiently accurate approximation for the practi- 
cal purpose of studying the electronic structure of solids. 
This may be understood from the fact that the same 
approximation, if applied to the uniform electron liquid, 
gives an excellent description of long-, intermediate-, and 
short-range correlations and furthermore obeys the Pauli 
principle; it can also give a highly accurate evaluation of 
the ground state energy for the electron liquid. This is 
because the self-energy functional employed is exact. 



We may say that Kohn-Sham equations assume the 
complete cancellation between the nonlocality (or the 
wavenumber dependence) and the energy dependence of 
£(r, r', e; [G]) in the evaluation of quasiparticle energies. 
This assumption is rather good, as we have already dis- 
cussed in Sec. V A. 

As for the energy gap of semiconductors, the gap of 
Kohn-Sham equations is related to that of Dyson equa- 
tions as E^ a p + A xc — E gap where A xc denotes the dis- 
cont i nuity of the exchange-correlation potential v xc (r) 
lllllllisf . According to an accurate Hartree-Fock calcu- 
lation in DFT HI, the relation E^ KS + A X = Eg? has 
been ascertained for silicon within an accuracy of 3%. It 
has been reported that the value of is a good 

approximation of the experimental gap and the disconti- 
nuity A x amounts to about twice the experimental gap. 
This implies that must be for the most part cancelled 
by A c . According to standard many-body perturbation 
theory, the well-known overestimation in the gap of semi- 
conductors calculated from usual Hartree-Fock equations 
can be remedied by the screened exchange potential. In 
the framework of DFT, on the other hand, the same 
screening effect can be represented by a great deal of 
cancellation between A x and A c . Probably, the total dis- 
continuity A xc (— A x + A c ) will be small but positive by 
analogy with the difference in the energy dispersion of the 
electron liquid between Dyson equations and Kohn-Sham 
equations. Kohn-Sham equations overestimate the den- 
sity of states at the Fermi level by less than 10% through 
the scaling factor m/m*. 



C. Correspondence between two self-consistencies 

In an attempt to interpret the extraordinarily en- 
hanced effective mass of heavy fermio n sys tems using the 
Hubbard model, a number of authors |ll3j | have exagger- 
ated the role of the energy dependence of the self-energy 
in the determination of the quasiparticle energy disper- 
sion around the Fermi level, neglecting its wavenumber 
dependence. In fact, there is a great deal of cancellation 
between the wavenumber dependence and the energy de- 
pendence of the self-energy. The former somewhat over- 
comes the latter in magnitude as we have elucidated in 
Sec. V A. This characteristic of the self-energy can be 
traced back to the dynamical nature of screening and 
hence is expected to apply to all metals as far as the 
Fermi surface is well defined. Then we have to ascribe 
the origin of the extraordinarily enhanced effective mass 
to a different charateristic of the self-energy in Dyson 
equations. 

In this subsection we would like to point out a char- 
acteristic of the self-energy peculiar to nonuniform real 
metals as a promising candidate for the origin of the ex- 
traordinarily enhanced effective mass. In transition met- 
als and rare earth metallic compounds the Fermi surface 
is constructed from a few kinds of electrons with differ- 
ent orbital characters, in contrast with the electron liquid 
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which consists of a single kind of electrons with the same 
character. 

Consider the self-energy S CT (r, r', e; [G]) entering Dyson 
equations. It is a nonlocal and energy-dependent opera- 
tor, but it is, at the same time, a functional of the self- 
consistent Green's function G. In the uniform electron 
liquid there is no essential difference in the neighborhood 
of the Fermi level between the first iterative and the fully 
self-consistent solutions for the quasiparticle energy dis- 
persion. In the case of rare earth metallic compounds, 
on the other hand, there is a very strong possibility that 
being a functional of the self-consistent Green's function 
G may cause a drastic change in the quasiparticle energy 
dispersion particularly in the immediate vicinity of the 
Fermi level due to the presence of f electrons and gives 
an extraordinarily sharp enhancement in the density of 
states. Probably, such a drastic change in the dispersion 
will be caused by the correlation-induced reconstruction 
of quasiparticle states in the immediate vicinity of the 
Fermi level and of course involve a significant change in 
the electron density n(r) of the system. 

According to Casida, the exchange-correlation poten- 
tial v xc (r) in Kohn-Sham equations can be interpreted to 
be the variationally best optimized local approximation 
to the self-energy E(r,r',e; [G]). Kohn-Sham equations 
assume the complete cancellation between the wavenum- 
ber dependence (non-locality) and the energy dependence 
of the self-energy, but still they maintain its essential 
property as a functional of the self-consistent G in the 
form of the exchange-correlation potential v xc (r) as a 
functional of the self-consistent electron density n(r). In 
other words, the implicit dependence of S(r, r', e; [G]) 
on the self-consistent G is in one-to-one correspondence 
to the implicit dependence of v xc (r; [n]) on the self- 
consistent n(r) in the sense that the latter occurs only 
from the former. This may be seen from the first iterative 
form of the Sham-Schlutcr equation which is identical to 
the OEP integral equation. 

/ dydy ' J ^ GKS ^y^) s (y-y') 

xv xc (y; [n])G K5 (y',r,c) 

= / J ^G^(r,y,c,)S(y,y>;[G^]) 

xG KS (y',r,c), (58) 

where G KS (r, y, uo) denotes the unperturbed Kohn-Sham 
Green's function in the many-body perturbation expan- 
sion. 

We may therefore expect that such a correlation ef- 
fect on quasiparticle energies as involves a change in 
the electron density n(r) of the system probably will 
be for the most part taken into account in the frame- 
work of Kohn-Sham equations in DFT. This is parallel 
to the expectation that there will probably be no signif- 
icant difference between the Kohn-Sham Fermi surface 
and the true Fermi surface, though strictly sp eaking the 
two Fermi surfaces are not identical |l!4l Ill5j | . The one- 



to-one correspondence between the two self-consistencies, 
i.e., S(r, r', e; [G^ s ]) as a functional of the self-consistent 
G KS and v xc (r; [n]) as a functional of the self-consistent 
n(r) will persist even for strongly anisotropic systems 
such as heavy fermion systems. 

In Sec. IV we have predicted that the correlation- 
induced reconstruction of Kohn-Sham orbitals in the 
vicinity of the Fermi level or the energy gap will be es- 
sential to the theoretical interpretation of the electronic 
structure of the so-called strongly correlated electron sys- 
tems. This prediction is based on the one-to-one corre- 
spondence between the two self-consistencies we have just 
mentioned above. 

To conclude, if the correlation energy functional to be 
used is sufficiently accurate, Kohn-Sham equations are 
expected to give a good approximation to quasiparticle 
energies as far as they do not exceed a critical energy 
near \x + hw v i. The orbital-dependent correlation energy 
functional we have proposed in Sec. IV B to be used 
in combination with the orbital-dependent exchange en- 
ergy functional probably will be sufficiently accurate to 
promote development in the band theory of solids. 



VI. CONCLUDING REMARKS 

One of the most important subjects in solid state 
physics at present is exploitation of more accurate meth- 
ods for treating correlation in the theoretical study of 
the electronic structure of solids. For this purpose we 
have constructed an orbital-dependent correlation energy 
functional E c [{ipi}] to be used in combination with the 
orbital-dependent exchange energy functional ^[{t/^}] 
in density-functional theory (DFT) with the help of de- 
tailed knowledge of long-, intermediate-, and short-range 
correlations of the electron liquid, not its density depen- 
dence of the correlation energy. This is a sophisticated 
but feasible method which may develop the band theory 
to a new advanced stage beyond the LDA. The present 
correlation energy functional £^[{1^}] in its construction 
is valid for tightly binding electrons as well as for nearly 
free electrons. 

Generally, effects of correlation on the electronic struc- 
ture of solids are so complicated and delicate that in order 
to obtain a reliable conclusion one has to take account 
of various aspects of correlation systematically starting 
from the original Coulomb interaction between valence 
electrons participating in cohesion of solids instead of re- 
sorting to intuitive simplifications such as the Hubbard 
model. 

Two aspects of long-range correlation arising from 
long-range parts of the Coulomb interaction in the elec- 
tron liquid, i.e., the existence of plasmons as a collective 
excitation and screening of the Coulomb interaction be- 
tween individual excitations (quasiparticles) as its com- 
plementary effect in the lower energy region apply to 
real metals though both are somewhat modified by band 
structure effects. The RPA is the simplest approxima- 
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tion to describe these aspects and its screening is in- 
dispensable to overcome the exchange difficulty. How- 
ever, it is important to recognize that the qualitative 
and quantitative evaluation of the quasiparticle energy 
dispersion E(p) and the effective mass m* in the region 
of metallic densities requires a systematic inclusion of 
higher-order perturbation terms obeying the Pauli prin- 
ciple in the expansion of the self-energy operator in the 
screened Coulomb interaction. This is because long-, 
intermediate-, and short-range correlations occurring in 
this region are very much complicated under the con- 
straint of the Pauli principle, in marked contrast to the 
situation in the high density region where long-range cor- 
relation in the RPA is predominant. In the evaluation 
of E(p) and m* of the electron liquid, a great deal of 
cancellation is observed between the wavenumber depen- 
dence and the energy dependence of the self-energy oper- 
ator. This can be interpreted as the dynamical nature of 
screening and hence we may expect that this cancellation 
will be observed in real metals so far as the Fermi surface 
is well defined. 

No qualitative difference may be expected in the neigh- 
borhood of the Fermi level or the energy gap between the 
quasiparticle energy dispersion calculated from Dyson 
equations and the one-electron energy dispersion from 
Kohn-Sham equations. The quantitative difference is less 
than 10% so far as we know. Then we may conclude that 
the use of Kohn-Sham equations will be valid for further 
development of the band theory. 

The reason for the unexpected success of the LDA may 
be ascribed to the long-range nature of the Coulomb in- 
teraction between electrons, which tends to cancel the 
difference between a real metal and its locally uniform 
electron liquid model in the quantitative evaluation of 
various electron-electron interaction effects. 

Except for the correlated properties coming from the 
electron liquid, what new correlation is the most im- 
portant in real metals at all? We answer this question 
from the viewpoint of DFT. The present -E c [{<^j}] is valid 
for both nearly-free and tightly-binding electrons since it 
assumes an orbital-dependent second-order perturbation 
like form in which all effects of higher-order perturbation 
terms are represented by an effective interaction defined 
for the electron liquid in the region of metallic densi- 
ties. A striking merit of the present i? c [{<pi}] is that 
it makes the meaning of correlation clear against the 
Hartree-Fock approximation in the framework of DFT 
and that it suggests how the electron density n(r) of the 
system may be redistributed under the influence of cor- 
relation. The present i? c [{</?i}] hybridizes Hartree-Fock 
Kohn-Sham orbitals in order to produce sclf-consistently 
reconstructed Kohn-Sham orbitals through the correla- 
tion potential v c (r). The important point is that the 
present v c (r) probably gives rise to less screening of the 
attractive nuclear potential at short distances to enhance 
the electron density n(r) in the immediate vicinity of in- 
dividual nuclei in solids. 

The electronic structure of the so-called strongly cor- 



related electron systems such as heavy fermions and 
Mott insulators might be interpreted in terms of sclf- 
consistently reconstructed Kohn-Sham orbitals in the 
neighborhood of the Fermi level or the energy gap. 

The expression "strongly correlated" has been exten- 
sively used in the sense that the electronic structure is 
beyond the scope of the LDA. In fact, the magnitude of 
the spin-antiparallel pair correlation function g^(r, r') 
for conduction electrons in cesium will be most reduced 
at short separations. Correctly, we should state that the 
LDA fails to give the proper description of correlation- 
induced reconstruction of Kohn-Sham orbitals. This 
reconstruction may be interpreted as a natural conse- 
quence of the fact that short-range correlation between 
valence electrons participating in cohesion of solids in- 
duces inhomogeneity in the electron density n(r) associ- 
ated with the presence of the nuclei, or equivalently short- 
range correlation between valence electrons significantly 
reduces the Hartrec- or Hartree-Fock screening of the at- 
tractive nuclear potential at short distances to enhance 
the electron density n(r) at the positions of individual 
nuclei in solids. This correlation-induced reconstruction 
of Kohn-Sham orbitals involving the enhancement of n(r) 
at the positions of individual nuclei is precisely the most 
important correlation that the LDA fails to describe. 

The present theoretical treatment of correlation will 
manifest the true value of DFT in the study of the 
electronic structure of solids. The present theory re- 
quires highly complicated self-consistency because of the 
orbital-dependence of the exchange-correlation energy 
functional E xc (= E x + E c ). Much technical difficulty is 
expected in its practical application. The most difficult 
problem will be the accuracy required in solving the OEP 
integral equation to evaluate the exchange-correlation 
potential v xc (r)(— v x (r) + v c (r)). The KLI method of 
approximately solving the OEP integral equation may be 
hopefully applied to the present orbital-dependent corre- 
lation energy functional. It is not so difficult to construct 
the spin-polarized version of the present theory. 
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